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SECTION  9  -  REFERENCES 


PREFACE 


by  M.  B.  McPherson 


ABOUT  THE  MODEL 

Without  doubt,  the  following  report  comprises  one  of  the  most 
significant  technical  memoranda  of  the  ASCE  Program.  Documented  is 
a  computer  model  that  should  see  extensive  use  in  total-jurisdiction 
preliminary  sewerage  planning,  if  for  no  reason  other  than  the  simple 
fact  that  it  is  presently  the  outstanding  tool  available  for  that  purpose. 
Combined  with  the  public  availability  of  the  user's  manual  (1)  and  the 
computer  program,  little  excuse  remains  for  any  local  government  agency 
not  to  investigate  the  potential  applicability  of  the  model  for  its  master 
planning. 


In  the  Preface  of  an  earlie  r  technical-memorandum  (2)  we  noted 
that  mathematical  models  used  for  the  simulation  of  urban  rainfall- runoff 
or  rainfall- runoff-quality  can  be  divided  into  three  distinct  categories: 
planning  models,  design  models  and  operations  models.  We  noted  that 
planning  models  are  used  in  massive  applications,  such  as  for  metro¬ 
politan  or  city- wide  master  plans.  As  an  example  of  the  scale  encountered 
the  City  of  Milwaukee  has  1,370-miles  of  separate  storm  drains  and 
combined  sewers  within  the  97-sq.  mi.  of  the  City  (3),  and  these  conduits 
are  distributed  over  465  drainage  catchments  having  a  maximum  size  of 
1,820-acres  and  a  median  size  of  25-acres  (4).  When  dealing  with  so  many 
components  the  model  used  must  be  as  simple  and  as  flexible  as  possible. 


(1)  Hydrologic  Engineering  Center,  Corps  of  Engineers,  "Urban  Storm 
Water  Runoff:  "STORM',"  Generalized  Computer  Program  723-S8- 
L2520,  Davis,  California,  May  1974, 

(2)  Lanyon,  Richard  F.  ,  and  James  P.  Jackson,  The  Metropolitan 
Sanitary  District  of  Greater  Chicago,  "A  Streamflow  Model  for 
Metropolitan  Planning  and  Design,"  ASCE  Urban  Water  Resources 
Research  Program,  Technical  Memorandum  No.  20,  ASCE,  New 
York,  N.  Y. ,  January  1974. 

(3)  Prawdzik,  Ted  B.  ,  Milwaukee  Department  of  Public  Works, 
"Environmental  and  Technical  Factors  for  Open  Drainage  Channels 
in  Milwaukee,  "  ASCE  Urban  Water  Resources  Research  Program, 
Technical  Memorandum  No.  12,  ASCE,  New  York,  N.Y., 

February  1970. 

(4)  Tucker,  L.  S.  ,  "Sewered  Drainage  Catchments  in  Major  Cities," 
ASCE  Urban  Water  Resources  Research  Program,  Technical 
Memorandum  No.  10,  ASCE,  New  York,  N.Y.,  March  31,  1969. 


That  is,  data  processing  for  planning  applications  becomes  a  much  more 
important  practical  consideration  than  the  level  of  sophistication  of 
hydrological  process  modeling,  whereas  just  the  opposite  emphasis  is 
required  for  design  applications. 

Design  models  can  be  very  elegant  and  detailed  tools  because 
they  are  used  for  analyzing  individual  catchments  and  subcatchments  in 
"one-shot"  applications  where  the  simulation  of  detailed  performance  of 
discrete  elements  within  a  subcatchment  must  be  achieved.  Whereas 
hourly  rainfall  data  is  an  appropriate  input  for  planning  models  and  for 
simulating  flows  in  larger  streams,  5-minute  interval  rainfall  data  (the 
shortest  duration  reported  by  the  U.  S.  Weather  Service)  is  the  appropriate 
input  for  simulating  flows  in  sewers  and  small  urban  streams  for  design 
applications.  Design  models  are  used  as  tactical  tools  and  planning  models 
are  used  as  tools  of  strategy;  and  operations  models  necessarily  embody 
both  of  these  capabilities. 


HISTORY  AND  ACKNOWLEDGMENT 

The  initial  version  of  the  model  was  employed  in  part  of  the 
development  of  the  Department  of  Public  Works,  City  and  County  of  San 
Francisco,  Master  Plan  for  combined  sewer  overflow  abatement  (5),  by 
Water  Resources  Engineers  (WRE).  In  February,  1973,  as  part  of  a 
training  course  on  "Management  of  Urban  Storm  Water,  Quantity  and 
Quality"  sponsored  by  The  Hydrologic  Engineering  Center  (HEC)  of  the 
Corps  of  Engineers,  an  advanced  version  of  the  model  was  presented  and 
about  one-third  of  the  total  course  time  was  devoted  to  hands-on  use 
of  the  model  by  participants.  WRE  developed  this  newer  version  and 
conducted  the  course  for  the  HEC.  Since  then,  the  HEC  developed  a 
user's  manual,  in  September  of  1973,  which  was  revised  in  May  of  1974 
(1).  Added  to  the  model  by  HEC,  and  included  in  the  explanations  of 
the  newer  manual,  are  capabilities  for  computing  quantity  and  quality  of 
rxmoff  from  nonurban  areas,  snowfall  and  snowmelt,  and  land  surface 
erosion  for  urban  and  nonurban  watersheds.  Because  explanations  of 
these  features  are  sufficiently  detailed  in  the  newer  manual  (1)  they  are 
not  repeated  in  the  following  report.  Dr.  Roesner  of  WRE  assembled 
the  following  report  with  help  from  the  coauthors,  using  in  part  the  texts 
of  three  HEC  course  'ectures  on  the  model.  (Eight  other  lectures  from 
the  course  comprise  a  companion  technical  memorandum  (6).  An  outline 


(5)  Department  of  Public  Works,  "San  Francisco  Master  Plan  for  Waste 
Water  Management,  "  City  and  County  of  San  Francisco,  September 
15,  1971,  (In  four  parts). 

(6)  Water  Resources  Engineers  and  the  Hydrologic  Engineering  Center  - 
Corps  of  Engineers,  "Management  of  Urban  Storm  Runoff,"  ASCE 
Urban  Water  Resources  Research  Program,  Technical  Memorandum 
No.  24,  ASCE,  New  York,  N.Y.,  May  1974. 


description  of  the  features  of  the  San  Francisco  Master  Plan  is  included 
in  a  previous  technical  memorandum  (7)*) 


We  are  greatly  indebted  to  the  authors  and  their  organizations 
for  the  privilege  of  presenting  this  unique,  break-through,  important 
report.  Indue  course,  copies  will  be  available  from  the  National  Technical 
Information  Service.  ASCE  Program  issuance  is  necessarily  restricted 
to  active  Program  cooperators,  numbering  over  two  hundred  persons. 
However,  not  only  are  Program  products  not  copyrighted,  but  anyone 
who  wishes  is  welcomed,  indeed  urged,  to  make  as  many  copies  as  they 
can  use  to  enhance  dissemination.  All  we  request  is  that  the  total  contents 
be  reproduced  to  insure  contextual  integrity. 

In  closing,  we  are  impelled  to  note  that  the  following  report 
epitomizes  a  central  Program  objective:  advancing  the  state  of  the  art. 


(7) 
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McPherson,  M.B.,  "Innovation:  A  Case  Study,"  ASCE  Urban  Water 
Resources  Research  Program,  Technical  Memorandum  No.  21, 
ASCE,  New  York,  N.  Y.  ,  February  1974. 


SECTION  1 
INTRODUCTION 


BACKGROUND 

It  has  only  been  within  the  last  decade  that  the  real  pollution 
potential  of  urban  runoff  has  come  to  be  recognized.  In  a  report  published 
in  1964  by  the  U.S.  Public  Health  Service  [l],  *  the  nationwide  signifi¬ 
cance  of  pollution  from  urban  runoff  was  first  identified.  Since  that  time, 
large  amounts  of  effort  and  money  have  been  devoted  to  the  characteri¬ 
zation  of  the  quality  of  urban  runoff  and  to  the  development  of 
methodologies  and  processes  to  control  this  source  of  pollution.  Funding 
for  these  studieshas  come  from  a  number  of  municipalities,  some  states, 
and  from  federal  agencies,  notably  the  Environmental  Protection  Agency 
(EPA)  and  recently  the  U.S.  Army  Corps  of  Engineers. 

A  review  of  these  studies  (see  References  10  and  11)  shows 
that  much  work  has  been  done  in  the  following  areas: 

1.  Development  of  stormwater  treatment  processes; 

2.  Sewer  system  control  to  maximize  pipeline  storage, 
thereby  reducing  the  amount  and  frequency  of 
overflows;  and 

3.  Characterization  of  the  quality  of  stormwater  and 
combined  sewer  discharges. 

In  addition,  several  sophisticated  mathematical  models  have  been 
developed  (some  with  funds  from  the  private  sector)  that  describe  the 
time-varying  hydraulic  response  of  an  urban  drainage  system  to  rainfall. 
A  few  of  these  models  include  descriptions  of  the  quality  of  urban  runoff. 
The  EPA  Stormwater  Management  Model  [2]  is  a  typical  example  of  the 
detail  and  scope  contained  within  these  models. 

The  information  and  technological  tools  that  are  presently 
available  are  deficient,  however,  in  that  they  do  not  adequately  address 
some  of  the  initial  questions  that  must  be  answered  in  the  preliminary 
planning  stage.  One  of  the  questions  that  must  be  answered  prior  to 
developing  a  pollution  control  plan  is:  what  is  the  present  and  expected 
future  magnitude  of  pollution  loads  carried  by  urban  runoff  from  a  given 
watershed?  Extensions  of  this  question  include  such  things  as  1)  what 
is  the  pollution  load  for  an  average  event,  2)  what  is  it  for  an  extreme 
event,  and  3)  how  often  does  a  given  extreme  event  occur? 


^'Numbers  in  brackets  refer  to  references  contained  in  Section  9. 


Given  that  we  can  answer  these  questions,  it  is  then  possible 
to  identify  some  constraints  for  the  stormwater  system  that  is  ultimately 
designed  so  that  the  receiving  waters  will  be  adequately  protected  (which 
is  the  whole  idea  of  "stormwater  management"  in  the  first  place). 
Therefore,  it  is  necessary  to  identify  those  systems  (i.e.,  combinations 
of  treatment  rate  and  storage  volume)  that  can  meet  the  constraints.  The 
extremes  of  these  combinations  are  obvious:  All  the  runoff  could  be  treated 
as  it  arrives  at  the  treatment  plant,  or  all  the  runoff  could  be  stored 
for  later  treatment  at  a  conventional  treatment  plant  during  off-peak 
hours.  Either  of  these  two  alternatives,  however,  will  normally  prove 
to  be  highly  uneconomical.  In  between  the  se  two  extremes  lie  an  extremely 
large  number  of  treatment  rates  and  storage  capacities  that  will  satisfy 
the  environmental  constraints  placed  on  the  system.  The  problem  is 
to  identify  the  feasible  combinations. 


REPORT  PURPOSE 

The  purpose  of  this  document  is  to  present  an  analytical  method 
that  can  be  used  in  the  preliminary  planning  stage  to  help  answer  the  types 
of  questions  posed  above.  The  method  has  been  coded  into  a  computer 
program  called  STORM  (Storage,  Treatment,  Overflow,  and  Runoff 
Model).  This  program  represents  a  method  of  analysis  to  estimate  the 
quantity  and  quality  of  runoff  from  small,  primarily  urban,  watersheds. 
Nonurban  areas  may  also  be  considered.  Land  surface  erosion  for  urban 
and  nonurban  areas  is  computed  in  addition  to  the  basic  water  quality 
parameters  of  suspended  and  settleable  solids,  biochemical  oxygen  demand 
(BOD),  total  nitrogen  (N),  and  orthophosphate  (P04).  The  purpose  of  the 
analysis  is  to  aid  in  the  selection  of  storage  and  treatment  facilities  to 
control  the  quantity  and  quality  of  urban  stormwater  runoff  and  land  surface 
erosion.  The  model  considers  the  interaction  of  eight  variables: 


1. 

precipitation, 

2. 

air  temperature  for  snowpack  accumulation 
and  snowmelt, 

3. 

runoff . "N 

4. 

pollutant  accumulation  1 

l  related  to 

on  the  land  surface .  i 

r  land  use 

5. 

land  surface  erosion,  .  .  .  J 

6. 

treatment  rates, 

7. 

storage,  and 

8. 

overflows  from  the  storage/treatment  system 

Land  uses  accounted  for  in  the  model  include:  single  family  residential, 
multiple  family  residential,  commercial,  industrial,  parks,  and  nonurban 
or  undeveloped  areas.  The  program  is  designed  for  use  with  many  years 
of  continuous  hourly  precipitation  records.  It  is  a  continuous  simulation 
model  but  may  be  used  for  selected  single  events. 

The  City  of  San  Francisco  used  this  program  in  the  preliminary 
planning  phase  of  their  Master  Plan  for  Stormwater  Management  [3].  The 
Corps  of  Engineers  is  currently  applying  STORM  in  several  of  their  Urban 
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Studies,  and  East  Bay  Municipal  Sewerage  District  No.  1  (which  serves 
seven  cities  on  the  east  side  of  San  F rancisco  Bay)  has  recently  used  the 
model  in  an  inflow/ infiltration  study  of  their  sanitary  sewer  system. 

The  program  has  been  documented  by  the  Hydrologic 
Engineering  Center,  Army  Corps  of  Engineers  in  Davis,  California,  and 
is  available  to  the  public  [4], 


HARDWARE  AND  SOFTWARE  REQUIREMENTS 

This  program  is  available  for  the  IBM  360/50,  UNIVAC  1108, 
and  CDC  6600  or  7600  computer  systems.  It  requires  about  35,  000  words 
of  core  storage  and  a  FORTRAN  IV  compiler  that  accepts  multiple  ENTRY 
statements.  Input  is  on  the  card  reader  and  possibly  a  tape/disk.  Output 
is  on  a  132  position  line  printer.  One  to  five  additional  tape/disk  units 
are  required  for  temporary  storage  during  the  processing.  The  only 
program  differences  among  the  three  computer  systems  are  due  to 
ENCODE/DECODE  type  statements  and  the  way  in  which  multiple  output 
files  are  handled.  Up  to  three  output  files  are  generated  on  tape/disk 
for  printing  at  the  end  of  the  job. 


REPORT  FORMAT 

The  concept  of  STORM  and  the  method  of  computing  runoff  is 
described  in  the  following  section.  In  Section  3  background  information 
on  the  quality  of  urban  runoff  is  presented  and  the  method  of  computing 
the  quality  of  urban  runoff  is  developed.  Section  4  discusses  computations 
of  treatment,  storage  and  overflow.  Section  5  describes  input  data 
requirements  for  STORM  and  the  output  information  it  produces.  Two 
actual  planning  applications  of  STORM  are  briefly  discussed  in  Section  6. 
Examples  of  other  possible  applications  of  STORM  are  contained  in  Section 
7.  Means  for  transferring  water  quality  data  for  one  location  to  another 
are  discussed  in  Section  8.  References  are  listed  in  Section  9. 
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SECTION  ? 

COMPUTATION  OF  RUNOFF  QUANTITY 


CONCEPT  OF  "STORM" 

The  quantity  of  urban  runoff  has  traditionally  been  estimated 
by  using  a  design  storm  through  frequency-duration-intensity  curves  or 
some  other  statistical  means  based  on  rainfall  records.  Such  approaches 
normally  neglect  the  spacing  between  storms  and  the  capacity  of  the 
urban  system  to  deal  with  some  types  of  storms  better  than  others. 

Often,  through  natural  and  artificial  storage  mechanisms, 
intense  short-duration  storms  maybe  completely  contained  within  storage 
so  that  no  untreated  stormwater  overflows  to  receiving  waters. 
Alternately,  a  series  of  closely  spaced,  moderately  sized  storms  may 
tax  the  system  to  the  point  that  excess  water  must  be  released  untreated. 
Consider,  for  example.  Figure  1  which  shows  the  response  of  two  different 
systems  to  the  same  ra  fall  trace.  System  A,  which  has  a  relatively 
high  treatment  rate  and  a  small  storage  capacity,  will  overflow  during 
the  high  intensity,  short  duration  storm.  However,  it  will  completely 
contain  the  second  storm  of  moderate  intensity  and  longer  duration.  System 
B,  on  the  other  hand,  which  has  a  low  treatment  rate  and  a  large  storage 
capacity,  completely  contains  the  first  storm.  Notice  that  it  would  also 
contain  the  second  storm  if  the  system  were  analyzed  independently  of 
the  antecedent  storm.  However,  in  this  case  the  spacing  of  the  storms 
is  such  that  the  system  analysis  must  include  both  rainstorms  as  a  single 
event  to  accurately  describe  the  system's  response  to  thd  rainfall  trace 
illustrated  in  the  figure. 

A  storm  cannot  be  defined  by  itself,  but  must  be  defined  taking 
into  account  the  response  characteristics  of  the  urban  stormwater  system 
It  is  for  this  reason  that  an  approach  was  developed  that  would  not  only 
recognize  the  properties  of  rainfall  duration  and  intensity,  but  would 
also  consider  storm  spacing  and  the  capacity  of  the  urban  stormwater 
system . 

Figure  2  shows,  pictorially,  the  interrelationship  of  the  eight 
stormwater  elements  considered  in  this  approach  for  estimating  storm¬ 
water  runoff  quality  and  quantity.  In  this  approach,  rainfall  washes  dust 
and  dirt  and  the  associated  pollutants  off  the  watershed  to  the  storage- 
treatment  facilities  so  that  as  much  stormwater  runoff  as  possible  can  be 
treated  prior  to  its  release.  Runoff  exceeding  the  capacity  of  the  treatment 
plant  is  stored  for  treatment  later.  When  the  storage  facilities  become 
inadequate  to  contain  the  runoff  the  untreated  excess  is  wasted  through 
overflow  directly  into  the  receiving  waters. 

For  a  given  precipitation  record,  the  quantity,  quality,  and 
number  of  overflows  will  vary  as  the  treatment  rate,  storage  capacity. 
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RAINFALL  TRACE 


System  Response  Exampl 


and  land  use  is  changed.  Land  surface  erosion  is  a  function  of  land 
use,  soil  types,  ground  slope,  rainfall/ snowmelt  energy  and  erosion 
control  practices.  A  typical  method  of  investigation  is  to  alter  the 
treatment,  storage,  and  land  use  and  note  the  resulting  response  of 
the  system.  A  group  of  alternatives  can  then  be  selected  from  among 
those  meeting  the  overflow  quantity  and  quality  objectives. 


COMPUTATION  OF  THE  QUANTITY  OF  RUNOFF 

Runoff  is  calculated  on  an  hourly  basis  as  a  function  of  rainfall 
plus  snowmelt  using  the  following  expression: 

R  =  C(P  -  f)  (1) 

where 

R  =  urban  area  runoff  in  inches  per  hour; 

C  =  composite  runoff  coefficient  dependent  on  urban 
land  use; 

P  -  rainfall  plus  snowmelt  in  inches  per  hour  over  the 
urban  area;  and 

f  -  available  urban  depression  storage  in  inches  per  hour. 

For  simplicity  we  will  omit  the  snowmelt  computation  in  our  discussion 
here.  The  interested  reader  is  referred  to  the  User's  Manual  [4]  for 
details  of  that  computation. 

The  runoff  coefficient  represents  losses  due  to  infiltration.  It 
is  computed  from  land  use  data  as  follows: 


C 


CP  +  (CI 


V 


m 

E 

L=  1 


XLFL 


(2) 


where 

Cp  =  runoff  coefficient  for  pervious  surfaces; 

Cj  =  runoff  coefficient  for  impervious  surfaces; 

X.  =  area  in  land  use  Las  a  fraction  of  total  watershed  area; 

Lj 

Fr  -  fraction  of  land  use  L  that  is  impervious;  and 
L 

m  =  total  number  of  urban  land  uses. 


.n 


A 


■  '  v-:-- 
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Before  the  runoff  coefficient  is  applied,  depression  storage 
losses  must  be  satisfied.  Depression  storage  represents  the  capacity  of 
the  watershed  to  retain  water  in  ditches,  depressions  and  on  foliage. 
The  amount  of  depression  storage  at  any  particular  time  is  a  function  of 
past  rainfall  plus  snowmelt  and  evapotranspiration  rates.  The  function 
is  computed  continuously  using  the  following  expression,  where  f  is  in 
inches: 

f  '  fo  +  N  D  k  *  for  f  -  D  O) 

where 

f  -  available  depression  storage,  in  inches,  after 
previous  rainfall; 

Np  -  number  of  dry  days  since  previous  rainfall; 

k  =  recession  factor,  in  inches/day,  representing  the 

recovery  (evapotranspiration)  of  depression  storage 
in  inches;  and 

D  -  maximum  available  depression  storage  in  inches. 

Figures  3a  and  3b  show  graphically  the  hourly  precipitation 
(P),  depression  storage  (f),  precipitation  excess  (P-f)  and  the  resulting 
runoff  (R).  Figures  3b  and  3c  show  how  the  runoff  is  distributed  between 
treatment  storage  and  overflow  for  a  system  with  a  treatment  rate  of 
0.02  inches/hour  and  a  storage  capacity  of  0.  16  inches. 


—  STORAGE  LEVEL 

—  STORAGE  CAPACITY 
■i  OVERFLOW 


FIGURE  3 

Time  Histories  of  Rainfall,  Runoff,  and  Storage 
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SECTION  3 

COMPUTATION  OF  RUNOFF  QUALITY 


SOURCES  OF  POLLUTANTS 

Basically,  pollutant  loads  are  introduced  into  urban  runoff  from 
three  sources: 

1.  The  land  surface  itself,  primarily  impervious  surfaces: 

2.  Catch  basins;  and 

3.  The  sewers  in  combined  systems. 

Of  these  three  sources,  the  land  is  the  most  important.  Catch  basins 
can  be  a  source  of  fi rst-flush  or  shock  pollution.  An  American  Public 
Works  Association  (APWA)  study  [5]  in  Chicago  found  that: 

.  the  liquid  remaining  in  a  basin  between  runoff 
events  tends  to  become  septic  and  that  the  solids  trapped 
in  the  basin  take  on  the  general  characteristics  of 
septic  or  anaerobic  sludge.  The  liquid  in  catch  basins 
is  displaced  by  fresh  runoff  water  in  the  ratio  of  one- 
half  the  volume  for  every  equal  volume  of  added  liquid. 

During  even  minor  rainfall  or  thaw  this  displacement 
factor  can  release  the  major  amount  of  the  retained 
liquid  and  some  solids.  The  catch  basin  liquid  was  found 
to  have  a  BOD  content  of  60  ppm  in  a  residential  area. 

For  even  minor  storms,  the  BOD  of  the  catch  basin  liquid 
would  be  seven-and-one-half  (7-1/2)  times  that  of  the 
runoff  which  had  been  in  contact  with  street  litter. 
Improved  design  of  catch  basins,  and  better  operational 
and  maintenance  practices,  could  reduce  this  first-flush 
pollutional  effect.  " 

In  combined  sewer  systems,  wastewater  is  incorporated  into 
the  storm  runoff.  In  addition,  the  storm  runoff,  as  it  passes  through 
large  sewers,  scours  sediment  deposited  by  wastewater  flows 
during  preceding  dry-weather  periods.  Figures  4  and  5  illustrate  the 
effects  of  wastewater  sewage  and  of  catch  basins  and  storm  sewer  scour 
on  the  quality  of  stormwater  overflows  [6], 

As  stated  above,  the  most  important  contributor  of  pollutants 
to  urban  runoff  is  the  land  surface  itself,  primarily  the  streets  and  gutters 
and  other  impervious  areas  directly  connected  to  streets  or  storm  sewers. 
Pollutants  accumulate  on  these  surfaces  in  a  variety  of  ways.  There  is. 
for  example:  debris  dropped  or  scattered  by  individuals:  sidewalk 

sweepings;  debris  and  pollutants  deposited  on  or  washed  into  streets  from 
yards  and  other  indigenous  open  areas;  wastes  and  dirt  from  building  and 
demolition;  fecal  droppings  from  dogs,  birds  and  other  animals;  remnants 
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FIGURE  5 

First-Flush  Pollutional  Effects  of  Catch  Basins 
And  Sewer  Scour  on  Combined  Sewer  Overflow 
(Selby  Street,  San  Francisco) 


of  household  refuse  dropped  during  collection  or  scattered  by  animals  or 
wind;  dirt,  oil,  tire  and  exhaust  residue  contributed  by  automobiles;  and 
fallout  of  air  pollution  particles.  The  list  could  go  on  and  on.  Irrespective 
of  the  way  in  which  pollutants  accumulate  on  the  urban  watershed,  they 
are  generally  associated  with  one  of  the  following  forms  of  street  litter: 

1 .  Rags , 

2.  Paper, 

3.  Dust  and  dirt, 

4.  Vegetation,  or 

5.  Inorganics. 

Table  1,  which  gives  estimated  street  litter  components  for  a  residential 
area  in  Chicago,  provides  a  rough  measure  of  the  relative  importance 
of  these  components. 


TABLE  1 

Monthly  Summary  of  Estimated  Street  Litter  Components, 
From  a  10-acre  (4  ha)  Residential  Area,  Chicago** 


Month 

Street  Refute  Component! 

(Ton*/ Month) 

Roy* 

Popei 

Dud  & 
Dirt 

Veqetntion 

Inoryonic 

Tn*ol  * 

Jon. 

.0015 

.OTA 
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.00 
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.0015 

.1  v» 
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,  v, 

.08 

.0V 

./*, 

Apr  I| 

.0016 

.ov. 

.65 

.Off 

.no 

.  //• 

Mo  y 

.  00 1 

.n.y, 

.56 

.00 

.00 

.  "r 

Jon# 

.ool  r. 

.ov, 

.ST 

.on 

.00 

. 

July 

.001  s 

,  OTA 

.56 

.on 

.0'' 

.  ft' 

Auy  , 

.001  s 

•  IU6 

..'>5 

.on 

.00 

.  v. 

5er»t. 

.001 r. 

.036 

.55 

,0H 

.07 

.  'A 

Oct. 

.oor. 

,0‘ftS 

.65 

.83 

.09 

*  .  '*/» 

Nov . 

.001  5 

.036 

.66 

.83 

.00 

1  .  >6 

Doc . 

.0015 

.036 

.55 

.00 

.09 

.68 

total* 

0180 

.43? 

6.60 

1.77 

1  .08 

10.4)1 

‘Some  total*  Hove  be<»n  rounded 


It  is  readily  apparent  that  the  most  significant  component  is  dust  and 
dirt  except  during  the  fall  of  the  year  when  vegetation  (primarily  leaves) 
becomes  the  dominant  component. 


""This  table  is  a  reproduction  of  Table  4  in  Reference  5. 


TYPES  OF  POLLUTANTS  AND  LOADING  RATES 


Nearly  all  of  the  pollutants  found  in  urban  runoff  are  associated 
with  the  dust  and  dirt  component  of  street  litter.  By  type,  COD,  BOD, 
and  solids  (suspended  and  settleables)  are  found  in  the  greatest  quantity. 
Nitrogen  and  phosphorus  are  also  found  in  significant  quantities.  In  areas 
where  street  deicing  by  salting  is  practiced,  winter  runoff  contains  very 
high  chloride  concentrations.  Other  pollutants  found  in  urban  runoff 
include  pesticides,  herbicides,  fertilizers  and  other  chemical  additives, 
heavy  metals,  and  many  other  known  and  unknown  pollutants. 

Data  on  the  rate  at  which  pollutants  accumulate  on  an  urban 
watershed  is  very  scanty.  In  fact,  it  is  almost  non-existent,  A  lot 
of  data  has  been  collected  on  the  quality  of  combined  sewer  overflows 
and  stormwater  discharges  for  various  cities  in  the  United  States  (Tulsa, 
Oklahoma;  Washington,  D.C.;  Atlanta,  Georgia;  San  Francisco, 
California;  Sacramento,  California;  and  Roanoke,  Virginia)  as  a  result 
of  the  U.S.  Environmental  Protection  Agency's  demonstration  grants 
program  for  abatement  of  stormwater  pollution.  The  studies  are  reported 
in  EPA' s  Water  Pollution  Research  Series;  however,  the  dissimilar  forms 
in  which  the  data  is  reported  makes  it  difficult  and  in  some  cases 
impossible  to  generalize.  The  problem  is  that  data  is  often  presented 
as  average  concentrations  or  as  pounds  of  pollutant  runoff  per  inch  of 
rain,  and  the  reported  values  may  be  for  combined  wastewater  and  storm 
runoff  rather  than  for  storm  runoff  alone. 

Even  for  a  given  watershed,  there  is  no  apt  description  of 
"typical"  stormwater  runoff  characteristics  because  of  the  variablility 
of  rainfall- runoff  patterns.  Thus,  reports  of  "mean  concentration"  or 
pounds  per  inch  of  rainfall  are  meaningless  as  generalized  variables  and 
they  show  poor  correlation  with  runoff  parameters. 

Results  from  a  demonstration  project  conducted  in  Tulsa, 
Oklahoma  [7],  were  summarized  in  terms  of  pounds  of  loading  per  day 
per  mile  of  street  for  each  of  15  areas  sampled  in  the  study.  These 
results,  presented  in  Table  Z,  give  an  indication  of  the  magnitude  of 
pollutant  buildup  for  different  land  uses.  These  findings  must  be  viewed 
with  caution,  however,  because  they  were  computed  by  taking  the 
"average"  concentration  of  the  pollutant  for  all  events  monitored,  which 
when  multiplied  by  the  total  storm  runoff  gave  "total  annual  mass 
emissions"  converted  into  a  rate  per  day  per  mile.  A  much  more  useful 
way  to  have  developed  this  data  would  have  been  to  sum  the  products 
of  discharge  and  concentrations  over  each  of  the  observed  events  and  then 
to  sum  all  such  total  storm  emissions  over  a  year. 

One  of  the  best  existing  sources  of  information  on  the  rate 
of  accumulation  of  pollutants  on  urban  watersheds  is  that  data  collected 
in  a  field  study  in  Chicago  by  APWA  [5],  This  study  determined  the  rate 
of  buildup  of  dust  and  dirt  in  the  streets  on  a  number  of  different  test 
areas  and  then  related  the  concentrations  of  various  pollutants  to  the  dust 
and  dirt.  A  summary  of  their  findings  is  contained  in  Table  3. 


TABLE  2 

Average  Daily  Loads  Per  Mile  of  Street 
(Tulsa,  Oklahoma) 


Total 

Test  Street 

Area  Miles 


Average  Load:  lbs/day/mile  of  street _ 

Organic 

Total  Kjeldahl  Soluble 

COD  Solids  Nitrogen  Orthophosphate 


Residential 

3 

14.  87 

1.41 

11.46 

120 

0.  26 

0.  34 

5 

16.  32 

2.  80 

21.43 

43 

0.  11 

0.  13 

7 

6.  84 

1.  20 

7.  20 

63 

0.  12 

0.  10 

8 

6.  97 

2.  72 

20.  89 

69 

0.  12 

0.  21 

9 

3.  11 

1.12 

13.  09 

47 

0.  07 

0.  11 

1 1 

49.  05 

1.  60 

13.  29 

66 

0.  08 

0.  15 

1  3 

5.  58 

2.  58 

15.  16 

81 

0.  25 

0.  20 

15 

2.  06 

2.47 

8.  67 

56 

0.07 

0.  17 

Comme  rical 

2 

7.  41 

2.  54 

15.  12 

92 

0.  32 

0.  29 

10 

12.  99 

2.  10 

20.  44 

82 

0.  16 

0.  13 

12 

3.  39 

4.  53 

25.47 

113 

0.  22 

0.  30 

Industrial 

1 

11.46 

4.  85 

41.  10 

838 

0.  41 

1.  30 

4 

28.  40 

3. 98 

29.  29 

175 

0/28 

0.  30 

6 

12.  24 

1 . 70 

12.  73 

49 

0.  09 

0.  1  3 

Average  Values 


Residential 
Commercial 
Industr  ial 


63.  1 
95.7 
354. 


Reproduced  from  Reference  7 


TABLE  3 

APWA  Findings  on  Rate  of  Pollutant  Buildup 
On  Urban  Watersheds* 


A'lmynt  of  Du»r  ond  o.f  ond  Strength  of  BOD  by 

Land  Use 
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L.iwl  IJ 

vt  1 1  •,  Jo f/  1  lK)  ft  of  curb 

_rEii _ 

■  ( i  *nm«* 

•  ..|>  3.3 

/./ 

1  lUUtl' 

.  4 .  6 

J 

*■  frtniify  .'..1 

3.6 

1 ,  <  et»deni.  e  0 .7 

') 

.•  'ed  oveiorie  1  > 

AitiouiiI  •_><  Poilv'ont  by  1  ype  of  Lo>>d  Use 

Item 

Single  family  Multiple  Family 
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.07 
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To  convert  the  data  ontained  in  Table  3  into  a  form  comparable 
with  that  of  Table  2,  the  Dust  and  Dirt  is  multiplied  by  2  (gutters  per 
street)  x  52.8  (100's  of  feet  of  gutter  per  mile)  x  (constituent  concen¬ 
tration/1000).  E.  g.  ,  the  rate  of  BOD  accumulation  on  an  urban  area  that 
is  single  family  residential  is:  0.7  x  (2  x  52.8)  x  5.0/1000  -  0.  3f 

Ibs/day/mile.  Rates  of  pollutant  buildup  in  pounds  per  day  per  mile  are 
given  in  Table  4. 


TABLE  4  . 

Average  Daily  Loads  Per  Mile  of  Street 
(Chicago,  Illinois) 


Average 

Load:  lbs 

/day/ mile 

of  street 

Land  Use 

BOD 

COD 

N 

P04 

Single  Family  Residential 

0.  36 

2.  95 

0.  03 

0.  004 

Multiple  Family  Residential 

0.  87 

9.  70 

0.  1  5 

0.012 

Commercial 

2. 70 

13.  6 

0.  14 

0.  024 

Industrial 

1. 45 

Comparison  of  the  values  in  Tables  2  and  4  indicates  that 
while  the  Chicago  data  is  consistently  lower  than  the  Tulsa  data,  the 
rates  of  buildup  are  similar  in  commercial  areas  except  for  P04.  Multiple 
family  residential  values  for  the  Chicago  area  also  compare  well  with 
the  Tulsa  data  except  for  P04.  The  principal  reason  for  the  lower  values 
in  the  Chicago  data  is  that  the  reported  values  are  for  the  soluble  portion 
of  the  constituents  and  do  not  include  that  portion  found  in  suspended  and 
settleable  material.  Also,  the  Tulsa  data  are  approximations  based 
on  over-simplified  computations.  Thus,  of  the  two  sets  of  reported  data 
the  Chicago  data  is  probably  better,  given  that  the  amount  of  constituents 
contained  in  the  solids  can  be  determined. 

One  fact  that  is  quite  evident  from  both  the  Chicago  and  Tulsa 
data  is  that  the  rate  of  buildup  of  pollutants  on  an  urban  watershed  varies 
significantly  with  land  use.  Intuition  would  tell  us  this  is  true.  Both 
sets  of  data  indicate  that  industrial  and  commercial  areas  are  much  dirtier 
than  residential  areas.  This  would  be  expected  since  there  is  higher 
pedestrian  and  vehicular  traffic  densities  in  these  areas.  The  data  shows 
that  pollutant  accumulation  rates  are  approximately  one  and  one-half  to 
five  times  as  great  in  commercial  and  industrial  areas  as  they  are  in 
residential  areas. 


ENTRY  OF  POLLUTANTS  INTO  URBAN  RUNOFF 

The  first  raindrops  that  fall  on  an  urban  watershed  simply  wet 
the  land  surface,  As  additional  rain  falls  the  impervious  surface  will 
become  wet  enough  that  some  of  the  water  begins  to  form  puddles,  filling 
the  depression  storage.  This  initial  rain  begins  to  dissolve  the  pollu¬ 
tants  in  the  gutters,  streets,  and  on  other  impervious  surfaces  and 
eventually,  as  this  water  actually  begins  to  flow  off  the  watershed  it 
carries  the  dissolved  material  in  it. 

As  rainfall  intensity  increases,  overland  flow  velocities 
become  sufficient  to  pick  up  solids.  Suspended  solids  are,  of  course, 
picked  upat  smaller  velocities  than  settleable  solids.  The  settleable  solids 
are  carried  off  the  watershed  in  two  ways.  If  the  velocity  is  sufficiently 
high,  the  settleable  solids  may  be  suspended  in  the  overland  flow.  At 
lower  velocities,  particles  may  simply  be  rolled  along  the  bottom  surface 
toward  the  stormwater  inlet. 

The  rain  that  initially  falls  on  pervious  surfaces  infiltrates 
into  the  ground.  If  the  rainfall  is  sufficiently  intense,  the  infiltration 
capacity  may  be  exceeded  and  the  excess  rainfall  begins  to  fill  the 
depression  storage  on  the  pervious  surfaces.  Finally,  if  the  rainfall 
is  of  sufficient  intensity  and  duration,  runoff  will  begin  to  flow  off  the 
pervious  areas,  onto  the  impervious  areas  and  thence  into  the  stormwater 
inlets.  Present  experience,  however,  indicates  that  the  amount  of  runoff, 
and  hence  the  pollution  loads,  contributed  from  pervious  surfaces  in  urban 
areas  are  small  compared  to  those  coming  from  the  impervious  areas 
and  can  be  neglected  in  determining  the  quality  of  surface  runoff.  This  is 
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especially  true  of  surfaces  covered  with  vegetation  such  as  lawns  and 
gardens.  Figure  6  illustrates  the  differences  in  runoff  and  pollution 
load  from  a  watershed  that  would  occur  if  it  was  converted  from  a  park 
(90%  pervious)  into  a  multiple  residential  area  (20%  pervious). 

ESTIMATION  OF  THE  RATE  OF  POLLUTANT  BUILDUP 
ON  URBAN  WATERSHEDS 

Since  dirt  is  the  major  component  of  street  litter  and  is  the 
primary  source  of  pollutants  in  urban  runoff,  the  most  basic  approach 
for  estimating  pollutant  buildup  rates  would  be  to  relate  them  to  the  dust 
and  dirt  accumulation  rates. 

Using  APWA  units  (Table  3)  for  the  rate  of  dust  and  dirt 
accumulation,  the  rate  of  buildup  DD^  for  a  given  land  use  L  can  be 
expressed  as: 

DDL  =  ddL  x  (G  L/100)  x  Al  (4) 

where 

DD^  =  rate  of  dust  and  dirt  accumulation  on  subareas  of 
land  use  L  in  lbs/day; 

dd  -  rate  of  dust  and  dirt  accumulation  for  land  use  L 
in  lbs/day/ 100  feet  of  gutter; 

G  =  feet  of  gutter  per  acre  for  land  use  L;  and 

L< 

A  j  =  area  in  land  use  L  in  acres. 

The  rate  factor  dd  should  be  supplied  by  the  user  for  his  area.  Default 

X-/ 

values,  which  are  those  shown  in  Table  3  are  incorporated  into  STORM 
and  can  be  used  if  no  better  data  are  available. 

The  initial  quantity  of  a  pollutant  p  on  subareas  of  land  use  L 
at  the  beginning  of  a  storm  can  then  be  computed  as: 

P  (F  x  DDr  x  N,_)  +  P  (5) 

p  p  L  D  po 

whe  re 

P  =  total  pounds  of  pollutant  p  on  land  use  L  at  the 
P  beginning  of  the  storm; 

F  -  pounds  of  pollutant  p  per  pound  of  dust  and  dirt; 
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N  j-j  =  number  of  dry  days  since  the  last  storm;  and 

P  =  total  pounds  of  pollutant  remaining  on  land  use  I, 
P°  and  the  end  of  the  last  storm. 


In  practice,  is  usually  limited  to  the  amount  that  would  be 

accumulated  in  a  90-day  dry  period.  The  reason  for  this  is  that  the 
efficacy  of  extrapolating  daily  buildup  rates  beyond  this  point  (which 
was  arbitrarily  selected)  is  uncertain.  Moreover,  if  equation  (5)  is  used 
repetitively  over  long  periods  of  time,  positive  errors  could  tend  to 
accumulate  in  P _ resulting  in  overly  large  values  of  P  . 


po 
If  street 


sweeping  is  practiced  on  the  watershed,  the  correct 


expression  to  use  for  P  is! 


=  P 


po 


(1  -E)  +  N  , 


DD, 


x  F  [( 1  -E) 
P  1 


( 1  -E)  n  1 


4  DD  r  X  F  (N  -  n  Nc  ) 
L  P  P  S' 


(6) 


whe  re 


Ng  =  number  of  days  between  street  sweepings; 

n  -  number  of  times  the  street  was  swept  since  the 
last  storm;  and 

E  =  efficiency  of  street  sweeping  (0,  6  to  0.  95). 


DETERMINATION  OF  URBAN  RUNOFF  POLLUTION  LOADS 


To  compute  the  amount  of  pollutant  washed  off  the  watershed 
during  a  storm,  it  is  assumed  that  the  amount  of  pollutant  removed  at 
any  time  t  is  proportional  to  the  amount  remaining: 


dP 

— E 


-KP 

P 


(7) 


We  stated  earlier  that  the  runoff  rate  Q  also  affects  rate  of  pollutant 
removal,  therefore  K  must  be  functionally  dependent  upon  Q.  However, 
given  two  identical  watersheds  except  for  their  area  size,  for  the  same 
rainfall  rate  r  on  both  watersheds  a  higher  runoff  rate  would  occur  from 
the  larger  watershed.  This  area  effect  can  be  eliminated  by  dividing  the 
runoff  Q  by  the  impe rvious  area  of  the  watershed.  The  impervious  area 
is  used  because  only  a  negligible  amount  of  the  runoff  comes  from  the 
pervious  area.  Since  cfs  per  acre  are  equivalent  to  inches  per  hour, 
we  can  say  that  K  is  functionally  dependent  on  the  runoff  rate  R  from 
the  impervious  area,  where  R  is  in  inches  per  hour.  Finally,  assuming 
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that  K  is  directly  proportional  to  R  and  that  a  uniform  rainfall  of  1/?. 
inch  per  hour  would  wash  away  90  percent  of  the  pollutant  in  one  hour 
(a  somewhat  arbitrary  assumption),  we  can  say  that  K  =  4.  6R.  Making 
this  substitution  into  equation  (7)  and  integrating  over  a  time  interval 
At  (during  which  R  is  held  constant)  gives: 

Pp  (t  +  At)  -  Pp  (t)e  '4‘  6RAt  (8) 

Equation  (8)  is  the  basic  form  of  the  overland  flow  quality  model  developed 
by  Metcalf  &  Eddy,  Inc.  ,  as  part  of  the  EPA  Stormwater  Management 
Model  [2],  Although  it  is  simplistic  and  contains  many  assumptions, 
it  is  the  best  overland  flow  water  quality  predictor  or  simulation  model 
that  presently  exists.  Moreover,  experience  with  that  model  (See 
Reference  6  and  Volume  II  of  Reference  2)  has  shown  it  to  give  fairly 
good  results. 

Some  idea  of  how  equation  (8)  behaves  can  be  gained  by  exami¬ 
nation  of  Figures  7  and  8.  Figure  7  shows  that  for  a  constant  runoff 
rate  R,  the  amount  of  pollutant  remaining  on  the  watershed  decays 
exponentially.  Under  a  time  varying  R  the  picture  is  quite  different,  as 
illustrated  by  the  upper  graph  Figure  8.  From  the  curve  of  P  vs.  t 
it  can  be  seen  that  the  amount  of  pollutant  removed  during  an  interval 
At  is  P(t)-P(t+  At).  The  rate  of  removal  of  mass  from  the  watershed  M 
is  simply  [P(t)-P(t+  A  t)  ]/  At,  which  can  be  expressed  as:  P 


Mp  =  P(t)  x  (1-e  '4l  6RA  S/At 


(9) 


The  variation  of  M  with  time 
P 


for  the  associated  hydrograph 


is  plotted  in  the  lower  graph  of  Figure  8.  A  plot  of  Mp  versus  t  is  termed 

a  pollutograph,  one  of  the  most  informative  methods  for  expressing  the 
pollutant  load  carried  by  urban  runoff.  To  determine  the  concentration 
of  a  pollutant  in  the  runoff  as  a  function  of  time,  one  simply  divides  the 
pollutograph  value  Mp  by  Q  (with  appropriate  conversion  factors). 

Equation  (9)  must  be  modified,  however,  because  not  all  of 
the  dust  and  dirt  on  the  watershed  is  available  for  inclusion  in  the  runoff 
at  a  given  time  t.  Thus  pollutants  which  are  tied  to  the  dust  and  dirt 
are  not  all  available  either.  The  Storm  Water  Management  Model  study 
[2]  found  that  for  suspended  solids  the  available  fraction  at  any  time  was: 


A  0.  0S7  +  1 . 4R 1  •  1  (10) 

SU8  '  1  u ' 

For  settleable  solids  it  has  been  assumed  that  the  availability  factor  is 

A  .  -  0.  028  +  1 .  OR  *’  8  (11) 

set  '  ' 

With  regard  to  BOD,  nitrogen  and  phosphate,  recall  that  the  APWA  data 
[S]  described  the  dissolved  fraction,  which  is  independent  of  the  amount  of 
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FIGURE  8 

Development  of  Pollutograph  (M  vs. 
F  rom  Time  History  of  P 


'■  V- 
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solids  available  for  runoff.  In  the  Storm  Water  Management  Model  study 
it  was  found  that  the  BOD  associated  with  the  suspended  solids  was  about 
10  percent  of  the  suspended  solids  load.  We  have  further  assumed  that 
the  BOD  tied  to  the  settleable  solids  is  2  percent  of  the  settleable  solids. 
For  nitrogen  and  phosphate,  we  have  assumed  that  BOD,  N,  and  P04 
are  associated  with  suspended  and  settleable  solids  in  the  same  proportion 
as  they  are  in  the  dissolved  state. 

Thus,  correcting  equation  (9)  for  available  suspended  and 
settleable  solids  and  adding  the  BOD,  N  and  P04  found  in  the  solids, 
we  get  the  following  set  of  equations  which  are  used  in  STORM: 


Suspended  Solids 


M  (t) 
sus 


=  A 


sus 


p  (t) 
sus 


EXPT 


where 


A  =  0.  057  +  1.4R1’  1 

SUS 

EXPT  =  ( 1 -e'4‘ 6RAt)/ At,  with  At 


1  hour 


(12) 


Settleable  Solids 


M  set(t> 


=  Aset  P.«t<« 


x  EXPT 


(13) 


where 

A  =  0. 028  +  R l* 8 
set 


BOD 


Mbod(t»  =  Pbod  «11  *  EXPT  +  10  M,us  4  °'02  U.U  "4I 

Nitrogen 

M  ..(t)  -  P  ..  (t)  x  EXPT  +  .045  M  +  .01  M  .  (15) 

nit  nit  sus  set 

P04 

MP04(t)  Fp04(t)xEXPT4  -0045Msus  +'001Mset<l6) 


SECTION  4 

COMPUTATION  OF  TREATMENT,  STORAGE  AND  OVERFLOW- 


PROCEDURE 

Computations  of  treatment,  storage  and  overflow  proceed  in 
an  hourly  step-by-step  method  throughout  a  period  of  rainfall / snowmelt 
record.  For  every  hour  in  which  runoff  occurs  the  treatment  facilities 
are  utilized  to  treat  as  much  runoff  as  possible.  When  the  runoff  rate 
exceeds  the  treatment  rate,  storage  is  utilized  to  contain  the  runoff. 
When  runoff  is  less  than  the  treatment  rate,  the  excess  treatment  rate 
is  utilized  to  diminish  the  storage  level.  If  the  storage  capacity  is 
exceeded,  all  excess  runoff  overflows  into  the  receiving  waters  and  does 
not  pass  through  the  storage  facility.  This  overflow  is  lost  from  the 
system  and  cannot  be  treated  later.  While  the  storm  runoff  is  in  storage 
its  age  is  increasing.  Various  methods  of  aging  are  used  including 
average,  first-in:  last-out,  first-in:  first-out,  or  others,  depending  on 
the  physical  conditions  encountered. 

The  computation  of  storage  and  the  interplays  among 
rainfall/snowmelt,  storage  and  treatment  represent  a  simplistic  approach 
for  dividing  a  rainfall  record  into  unique  events  such  that  the  event  is 
defined  in  terms  of  the  urban  system.  For  example,  whether  two  "storms" 
are  considered  as  two  isolated  occurrences  or  as  one  large  storm  is 
entirely  dependent  upon  how  the  system  will  react  to  them.  If  the  system 
has  not  recoverea  from  the  fi rst  when  the  second  arrive s,  the  two  definitely 
will  interact  and  hence  must  be  considered  together.  "Events"  are  defined 
as  beginning  when  storage  is  required  and  continues  until  the  storage 
reservoir  is  emptied.  All  the  rainfall  occurring  within  this  period  is 
regarded  as  oart  of  the  same  event.  If  precipitation  produces  runoff  that 
does  not  exceed  the  treatment  rate,  the  runoff  will  pass  through  the 
treatment  process  but  will  not  register  as  an  event.  From  the  standpoint 
of  the  urban  stormwater  system,  such  pre cipitation  is  inconsequential  and 
hence  is  not  part  of  an  "event"  even  if  it  should  occur  immediately  preceding 
an  obvious  event. 

The  runoff  coming  into  the  storage/treatment  system  is  given 
by  equation  (1).  The  quantity  of  system  overflows  are  computed  using: 


Q 

O 

=  R  -  Qt  -  Qs 

(17a) 

qt 

-  minimum  of  (R  +  Q  ,  T) 

St-1 

(17b) 

Qs 

-  minimum  of  (R  -  ,  S) 

(17c) 
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where 


Q 


s 


t-1 


watershed 

watershed 

watershed 

watershed 


runoff  overflow,  in  inches; 

runoff  treated,  in  inches; 

runoff  stored,  in  inches; 

storage  remaining  in  previous  hour, 


inches; 


R  =  watershed  runoff  as  calculated  using  equation  (1),  inches; 

T  =  treatment  rate  in  watershed  inches/hour;  and 
S  =  storage  capacity  in  watershed,  inches. 

The  quality  of  system  overflows  are  computed  as  follows  for  each  pollutant 

(18) 
(19) 


for  each  hour: 

M 

=  M  (Q  /R) 

po 

P  ° 

MpS/T 

=  M  -  M 

P  Po 

where 


M 


po 


total  pounds  of  pollutant  overflowing  from  system; 


M 


M 


P 

pS/T 


total  pounds  of  pollutant  p  coming  into  the  system;  and 
total  pounds  of  pollutant  p  going  to  storage/treatment. 


The  program  does  not  model  the  treatment  process  but  it  does 
compute  the  quantity  of  water  treated.  It  is  assumed  that  the  pollutants 
will  be  reduced  to  an  acceptable  level  before  the  storm  water  is  released. 
The  age  of  pollutant  in  storage  is  computed  as  previously  mentioned. 
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SECTION  5 

INPUT  AND  OUTPUT  FOR  "STORM" 


The  basic  input  data  required  by  STORM  and  the  basic  output 
data  generated  by  the  program  are  illustrated  in  Figure  9.  This  section 
defines  the  input  data  requirements  more  specifically  and  contains  details 
on  STORM  output. 


INPUT  DATA  REQUIREMENTS 
Hydrogeometric  Data 

The  first  step  in  setting  up  data  for  the  simulation  model  is  to 
define  the  boundaries  of  the  basin  which  is  to  be  investigated,  specifically 
that  area  which  drains  to  some  specific  point  of  inte  rest  such  as  a  receiving 
water.  The  size  of  the  area  is  a  computation  variable  but  it  should  be 
limited  to  less  than  10  square  miles  so  that  travel  time  in  the  system 
can  be  neglected. 

Once  the  drainage  basin  boundaries  are  set  the  following 
information  is  required: 

1.  Size  of  the  total  area  of  the  basin 

2.  Percent  of  the  total  area  in  each  of  the  following  land 
use  groups: 

a.  Single  Family  Residential 

b.  Multiple  Family  Residential 

c.  Commercial 

d.  Industrial 

e.  Open  or  park 

3.  Average  percent  imperviousness  of  each  land 
use  group 

4.  Feet  of  gutter  per  acre  for  each  land  use  group 

*>.  A  runoff  coefficient  for  impervious  areas  (the 
usual  range  is  0.  8  to  0.  9) 

f).  A  runoff  coefficient  for  pervious  areas  (the  usual 
range  is  0.  1  to  0.  3) 

7.  The  depression  storage  available  on  the  impervious 
areas  (usually  0.05  to  0.  1  inches). 
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Determination  of  the  percent  of  area  under  the  various  land 
uses  can  be  a  tedious  task.  However,  most  jurisdictions  have  this 
information  already  available  in  one  form  or  another,  or  they  have  maps 
of  sufficient  scale  that  the  various  land  uses  can  be  identified  and  their 
areas  calculated. 

Determination  of  the  average  percent  imperviousness  of  each 
land  use  group  must  be  done  carefully,  for  this  is  the  most  sensitive 
parameter  affecting  the  amount  of  storm  runoff  which  comes  off  a 
watershed.  In  particular,  impervious  surface  areas  which  drain  to 
pervious  areas  should  be  excluded  from  the  impervious  fraction  because 
runoff  from  these  surfaces  will  probably  be  held  on  the  pervious  surfaces 
to  which  they  drain.  In  this  regard,  special  attention  should  be  given  to 
whether  the  house  gutter  drains  to  pervious  land  or  to  the  sewer  system 
or  gutter.  Also,  in  residential  areas  where  there  are  parking  strips, 
the  sidewalks  will  most  likely  drain  to  the  pervious  area  on  either  side 
of  the  walk.  Table  5  may  offer  some  guidance  for  determining  the  percent 
imperviousness  for  different  land  uses. 


TABLE  5 

Percent  Imperviousness  for 
Various  Land  Uses  in  the  San  Francisco  Area* 


Type  of  development 

(1) 

Density, 
in  units 
per  acre 
(2) 

Percent  impervious 

Santa  Clara 
County 
(3) 

San  Francisco 
Bay  Region 
(4) 

Residential: 

Hill  areas 

0.5-  2 

6 

8 

Low  urbanization 

3-6 

10 

15 

Medium  urbanization 

7  -10 

20 

25 

Heavy  urbanization 

11  -20 

32 

40 

(apartments) 

Industr lal : 

Nonmanufacturing 

50 

60 

Manufacturing 

40 

50 

Reserve 

20 

25 

Corancrc i a  1 

50 

60 

Transpora  t ion 

70 

75 

Public  buildings 

40 

50 

Public  parks 

12 

12 

Agrlcu 1 tural 

4 

4 

Natural  watersheds 

2 

2 

:  Extracted  from  Table  1  of  Reference  8. 


Estimation  of  the  number  of  feet  of  gutter  per  acre  is  best 
done  from  a  plat.  But  if  none  is  available  a  reasonable  estimate  can  be 
based  on  the  average  size  of  the  block  by  taking  the  perimeter  of  the  block 
times  the  number  of  blocks  per  acre  from  a  street  map. 

Hydrologic  Data 


A  record  of  hourly  rainfall  is  required.  The  rainfall  record 


may  be  as  long  or  as  short  as  desired  but  should  be  of  sufficient  length 
to  assure  that  all  storms  of  interest  are  included  in  the  record.  Ten 
to  thirty  years  of  record  is  desirable.  A  long  raingage  record  exists 
for  most  cities.  Where  such  information  is  lacking,  however,  standard 
hydrologic  procedures  for  areal  translation  of  rainfall  records  will  have 
to  be  applied. 

Quality  Data 

The 

quality  data  required  for  the  simulation  model  consists  of: 

1. 

The  daily  rate  of  dust  and  dirt  accumulation  in 
pounds  per  100  feet  of  gutter  for  each  of  the  land 
use  areas: 

a.  Single  Family  Residential 

b.  Multiple  Family  Residential 

c.  Commercial 

d.  Industrial 

e.  Open  or  Park 

2. 

The  pounds  of  each  of  the  following  pollutants  per 

100  pounds  of  dust  and  dirt  for  each  land  use  category: 

a.  Suspended  solids 

b.  Settleable  solids 

c.  Soluble  BOD 

d.  Soluble  N 

e.  Soluble  P04 

3. 

The  interval  in  days  between  street  sweepings  for 
each  land  use  category 

4. 

Street  sweeping  efficiency  (usual  range  is  .  6  to  .  9). 

Because  this  data  is  difficult  to  obtain,  default  values  are 
provided  in  the  computer  program  as  follows: 
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1.  Daily  Rate  of  Dust  and  Dirt  Accumulation* 

Land  Use  Amount  of  D/D  by  Land  Use, 

lb./day/lOO  ft.  of  Gutter 


Single  Family  Residential  0.7 

Multiple  Family  Residential  2.3 

Commercial  3.  3 

Industrial  4.  f> 

Open  or  Park*5'-  1.5 


»  2.  Pounds  of  Pollutant  in  Dust 

r 

and  Di  rt* 

Land  Use 

Lbs. 

of  Pollutant/  100  lbs 

.  of  D/D 

Sus , 

Sett. 

Solid  s* 

*  Solids** 

BOD 

N 

P04  j 

Single  Family  Residential 

11.1 

1  .  1 

0.  5 

0.048 

0.005  I 

1  Multiple  Family  Residential  8.0 

0.  8 

0.  36 

0.  061 

0.005  j 

Commercial 

17.0 

1.7 

0.  77 

0.041 

0.007 

|  Industrial** 

6.  7 

0.  7 

0.  3 

0.04  3 

0.003 

1  Open  or  Park** 

11.  1 

1 .  1 

0.  5 

0.  048 

0.005  ! 

!  3.  Street  Sweeping 

Inte rval 

90  days 

i 

\ 

I  4.  Street  Sweeping 

Efficiency 

0.  7 

! 

!  1 

OUTPUT  FROM  "STORM" 

The  computer  program  produces  four  output  reports: 

1.  Quantity  Analysis, 

2.  Quality  Analysis, 

3.  Pollutograph  Analysis,  and 

4.  Land  Surface  Erosion  Analysis, 

For  the  quantity  and  quality  analyses,  STORM  gene  rates  statistics  by  event 
plus  the  average  statistics  for  all  events.  A  complete  list  of  the  output 
statistics  from  the  quantity  and  quality  analyses  are  contained  in  Table  6, 

Tables  7  through  11  are  examples  of  STORM  output  from  the 
Quantity  Analysis,  Quality  Analysis  and  Pollutograph  Analysis.  For 
details  on  the  output  from  the  Land  Surface  Erosion  Analysis  the  reader 
is  referred  to  the  program  manual  [4], 


"Data  is  taken  from  APWA  Chicago  Study  (See  Reference  5). 
"’"Estimated  values  from  other  sources. 


I 
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TABLE  (> 
"STORM”  Output 


I.  STATISTICS  BY  EVENTS 

A.  RAINFALL 

/  DURATION  OF  RAINFALL  EVENT 
2  HOURS  OF  RAIN 
3.  TOTAL  RAINFALL 

B.  STORAGE 

/  TIME  SINCE  LAST  EVENT 

2  DURATION  OF  STORAGE 

3  TIME  TO  EMPTY 

4  MAXIMUM  STORAGE  USED 

C.  OVERFLOW 

/  TIME  OVERFLOW  STARTS 

2  DURATION  OF  OVERFLOW 

3  quantity  of  ov,  clow 

4  OVERFLOW  IN  FIRST  THREE  HOURS 

D.  TREATMENT 

/  DURATION  OF  TREATMENT 

2.  QUANTITY  TREATED 

E.  QUALITY  (susp  solids,  sett  solids,  BOD,  nitrogen,  phosphorous) 
/  MASS  EMISSION  IN  RUNOFF 

2  MASS  EMISSION  OF  OVERFLOW 

3.  MASS  EMISSION  DURING 

FIRST  THREE  HOURS  OF  OVERFLOW 

I.  AVERAGE  STATISTICS  ( A- E  ABOVE) 

A.  FOR  ALL  EVENTS 

B.  FOR  ALL  OVERFLOW  EVENTS 

C.  EVENTS  /  YR 

D.  OVERFLOWS  /  YR 
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F*CE  t  Sf  L8Y  STREET  NaTERSHEO  SAN  FRANCISCO,  CA, 

.  quantity  analysis 

TREATMENT  rate  •  IN/1R,  37,8  CIS,  24.400  NGO  FED  BLOC,  SaN  FRANCISCO,  CA 

STORAGE  CAPACITY.  .1UJ3  INCNES,  31.2  AC-FT,  10.168  RG  SELBY  STREET 
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Example  of  Quantity  Statistics  by  E 


_  ,  TABLE  8 

Example  of  Average  Quantity  Statistics  for  all  Events 
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Example  of  Quality  Statistics  by  Event 


v  ,  TABLE  10 

Example  of  Average  Quality  Statistics  for  all  Events 
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SECTION  6 

PLANNING  APPLICATIONS 


CITY  OF  SAN  FRANCISCO 

The  City  of  San  F  rancisco  funded  the  original  development  of 
the  Quantity  Analysis  portion  of  STORM.  The  purpose  of  the  program  was 
to  create  a  tool  that  would  enable  the  City  to  evaluate  the  effectiveness 
of  various  combinations  of  treatment  rate  and  storage  capacity  with 
respect  to  their  ability  to  reduce  combined  sewer  system  overflows. 
Results  from  STORM  were  used  as  a  guide  in  the  initial  sizing  of  facilities 
for  the  City's  Master  Plan  [3],  The  following  paragraphs  are  abstracted 
from  the  Master  Plan  Report.  Additional  details  are  described  in 
References  12  and  13. 

The  rainfall  record  used  for  the  analysis  was  a  6?  year  U.S. 
Weather  Service  record  of  hourly  values  of  rainfall  measured  at  the 
Federal  Office  Building  in  the  City.  A  runoff  coefficient  of  0.65  was 
assumed  for  the  analysis.  Table  12  shows  the  various  combinations 
of  treatment  rate  and  storage  capacity  that  were  examined,  the  resulting 
events  per  year,  overflows  per  year  and  average  quantity  of  overflow  per 
year.  These  data  are  displayed  graphically  in  Figures  10  and  11,  which 
show  the  relationshipbetween  given  combinations  of  storage  and  treatment 
with  overflow  frequencies  and  with  overflow  volumes,  respectively. 

Application  of  STORM  to  the  Federal  Office  Building  record 
with  0  storage  and  a  treatment  rate  of  0.02  inches  per  hour  provided 
the  baseline  or  existing  condition  data.  From  this  computation  it  was 
determined  that  approximately  one -third  of  the  runoff  is  presently  treated 
and  discharged  by  the  three  water  pollution  control  plants  and  that  the 
other  two-thirds,  or  about  6.0  billion  gallons  of  runoff  per  year,  overflows 
without  treatment.  This  volume  of  overflow  occurs  during  an  aggregate 
average  of  206  hours  per  year.  On  the  average,  there  are  46  days  in  the 
year  during  which  8?  overflows  occur. 

The  storage  needed  to  contain  all  overflows  from  the  greatest 
recorded  storm  utilizing  the  existing  treatment  rates  would  be  240  million 
cubic  feet.  This  storage  volume  is  then  the  upper  limit  of  an  all-storage 
scheme  and  exceeds  by  a  factor  of  2  the  volume  requirement  of  an  all- 
treatment  scheme. 

The  data  presented  above  was  necessarily  based  on  two 
assumptions:  that  the  runoff  loss  is  35  percent  and  that  rainfall 
occurrence  is  uniform  over  the  City.  Each  is  a  significant  parameter 
in  determining  the  total  volumes  of  runoff.  At  the  time  the  Master  Plan 
was  being  developed  no  verified  data  existed  on  the  losses  experienced 
in  the  rainfall  -  runoff  process,  although  some  measurements  have  been 
made  in  more  recent  characterization  studies. 
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FIGURE  10 
Overflow  Frequencies 
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STORAGE  CAPACITY  -  INCHES 


FIGURE  1  1 
Overflow  Volumes 
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Storms  monitored  during  the  1969/1970  rainy  season  by  a 
system  of  19  rain  gages  distributed  over  the  city  showed  a  15  percent 
lower  overall  average  volume  of  rainfall  over  the  whole  city  than  that 
indicated  by  the  gage  at  the  Federal  Office  Building.  Because  the  time 
correlation  for  the  19  rain  gages  which  were  operational  during  that 
season  was  poor,  the  above  percentage  indicates  only  the  extent  of  spatial 
variation. 


Spatial  and  temporal  differences  observed  in  the  occurrence 
of  rainfall  led  to  the  conclusion  that  a  system  featuring  interconnection 
would  result  in  more  efficient  utilization  of  facilities.  Further,  the  use 
of  real-time  computer-actuated  control,  based  on  sensing  storm  direction 
and  the  likely  volumes  of  rainfall,  with  a  constant  concurrent  updating 
of  the  status  of  the  system,  would  permit  a  maximum  use  of  all  capacity 
throughout  the  system.  The  result  would  be  the  construction  of  fewer 
and  smaller  facilities  which  would  serve  the  overall  system  rather 
than  only  discrete  segments  of  the  system. 

Other  STORM  statistics  for  the  various  combinations  of  storage 
and  treatment  include  the  number  of  events  that  would  have  occurred,  the 
volume  of  overflows,  the  duration  of  overflows,  and  the  number  of  days 
of  overflows.  With  the  previous  assumptions,  this  data  was  used  for 
preliminary  evaluation  of  proposed  control  systems  with  regard  to  possible 
overflow  quality  and  mass  emissions  of  constituents. 

The  initial  sizing  of  the  Master  Plan  system  was  based  upon 
the  records  available  from  the  Weather  Service  gage.  This  data 
represented  the  best  information  available  at  the  time  and  all  other  data 
indicated  that  any  design  based  upon  this  gage  would  likely  be  conservative 
with  regard  to  size  and  costs.  Refinement  of  the  design  will  take  place 
as  data  accumulates  from  the  City's  extensive  field  information  collection 
system . 


Figure  12  shows  a  composite  of  the  effects  of  various 
combinations  of  storage  and  treatment  with  regard  to  the  frequency  of 
uncontrolled  overflow  occurrences.  It  is  apparent  that,  given  a  desired 
frequency  of  overflow  occurrence,  increasing  the  treatment  rate  decreases 
the  storage  requirements  and  that,  for  any  gi-ven  storage  volume, 
increasing  treatment  capacity  results  in  a  lower  occurrence  frequency. 
Further  inspection  of  the  figure  also  indicates  that  the  law  of  diminishing 
returns  results  in  increasingly  greater  storage  requirements  for  any 
treatment  rate  to  attain  a  lowe  r  f requency  of  ove r flow  occurrence .  Through 
the  application  of  i  ost  factors  for  storage  and  treatment  facilities, 
optimum  design  points  for  minimum  cost  for  various  levels  of  control 
can  be  derived. 

Sewe  r  Sy  stem 

The  present  storm  sewer  design  criteria  includes  the 
conveyance  of  a  5-year  intensity  rainfall  without  flooding.  When  rainfall 
intensity  exceeds  the  design  rate,  surface  transport  and  flooding  can 
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NUMBER  OF  OVERFLOW  /  YEAR 


occur.  As  may  be  seen  on  Figure  1"?,  there  are  numerous  locations  in 
the  city  where  surface  waters  can  accumulate  until  capacity  in  the  system 
will  accept  them.  If  detention  basins  were  located  at  these  sites  with 
appropriate  street  drains,  surface  waters  would  flow  to  the  basins.  Such 
facilities  installed  in  these  locations  should  provide  for  greater  public 
protection  from  the  inconvenience  of  surface  ponding. 

Another  benefit  of  detention  basins  to  the  conveyance  system 
is  obtained  by  limiting  the  flow  to  downstream  conduits. 

The  location  of  stormwater  detention  or  storage  basins  in  any 
particular  sewer  system  has  a  beneficial  effect  on  the  available  downstream 
conduit  transport  capacity  in  terms  of  historical  rates  of  rainfall.  The 
selection  of  storage  basin  locations  within  any  stormwater  conveyance 
system  can  be  made  such  that  the  main  trunk  sewers  of  the  system  down¬ 
stream  of  the  basins  will  be  upgraded  with  regard  to  the  size  of  the  rainfall 
that  may  be  conveyed  before  exceeding  sewer  capacities.  The  volume  of 
storage  facilities  considered  must  satisfy  the  following  criteria: 

1.  The  storage  volumes  utilized  shall  not  be  less  than 
1.2  million  cubic  feet.  This  restriction  stems 
from  the  consideration  of  economics  of  con¬ 
struction  of  such  basins. 

2.  The  required  volume  for  storage  is  equivalent  to 
the  volume  of  runoff  derived  from  1  inch  of  rain¬ 
fall  on  the  contributing  fraction  of  the  watershed 
which  is  not  already  tributary  to  the  basin.  This 
is  based  upon  the  frequency  of  overflow  of  such 
basins  which  historically  would  occur  when  used 
in  conjunction  with  the  next  criterion, 

3.  The  evacuation  of  flow  out  of  storage  is  continuous 
and  is  equivalent  to  the  rate  of  runoff  from  a 
steady  state  rate  of  rainfall  of  0,  10  in/hr  from 
all  upstream  tributary  areas.  This  criterion,  in 
conjunction  with  2  above,  provides  a  storage 
volume  which  historically  will  overflow  only  one 
time  every  five  years. 

It  is  estimated  that  the  cost  would  be  $77  million  to  replace 
all  sewers  of  inadequate  capacity  that  are  larger  than  three  feet  in  size. 
The  cost  of  basins  will  be  offset  to  some  degree  by  the  equivalent  costs 
foregone  for  installing  additional  or  longer  sewers  in  those  areas  where 
inadequacies  now  exist.  All  conceptual  design  costs  were  evaluated  with 
regard  to  this  aspect.  Of  $150  million  required  to  remove  the  inadequacy 
of  all  sewers,  about  $50  million  could  be  saved  by  means  of  the  detention 
system. 
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A  further  benefit  of  detention  basins  is  in  the  potential  for 
flushing  the  conveyance  system  with  storage  flows.  This  may  reduce 
maintenance  costs  for  lower  portions  of  the  system  that  are  in 
subsidence  areas. 


EAST  BAY  MUNICIPAL  UTILITY  DISTRICT 

The  East  Bay  Municipal  Utility  District,  Special  District  No.  1 
(EBMUD,  SD1),  provides  wastewater  disposal  services  to  seven  cities 
on  the  east  side  of  San  Francisco  Bay,  California.  Figure  14  shows  the 
location  of  the  service  area.  The  EBMUD  system  is  comprised  of  a  22  mile 
interceptor  system  plus  a  central  treatment  plant. 

Wet  weather  inflow/infiltration  is  a  major  problem  for 
EBMUD.  The  source  of  the  problem  is  the  wastewater  sewer  systems  that 
drain  into  the  interceptor  system.  Most  of  these  systerrsare  over  30  years 
old.  A  number  of  them  have  been  converted  from  combined  sewers  to 
"sanitary"  sewers,  and  a  few  combined  sewers  are  still  connected  to  the 
interceptor  system.  In  addition,  there  are  many  illegal  connections 
(especially  to  roof  and  yard  drains).  Present  estimates  are  that  11  percent 
of  the  rainfall  on  the  service  area  appears  as  inflow/infiltration  in  the 
interceptor  system  causing  the  system  to  overflow  to  San  Francisco  Bay 
about  11  times  a  year. 

Reduction  of  the  overflows  could  be  accomplished  by  reducing 
extraneous  inflows,  by  providing  additional  treatment  and  storage,  or 
through  some  combination  of  both.  The  question  was,  which  particular 
combination  would  give  the  control  required  and  which  combination  would 
be  the  most  cost-effective. 

It  was  estimated  by  the  District  staff  that  conversion  of  the 
remaining  combined  sewers  to  sanitary  sewers  would  reduce  the  gross 
inflow/ infiltration  ratio  (i.e.,  runoff  coefficient)  from  11  percent  to  8 
percent.  If  80  percent  of  the  direct  connections  to  the  sewer  system 
could  be  eliminated  (roof  and  yard  drains,  parking  lots,  catch  basins 
connected  in  error,  etc.  )  it  was  estimated  that  the  inflow/infiltration  ratio 
could  be  reduced  another  3  percent,  i.e.  from  8  percent  to  5  percent. 
Finally,  by  the  additional  removal  of  fifty  percent  of  the  percolation 
infiltration,  it  was  estimated  that  the  inflow/infiltration  ratio  could  be 
reduced  another  2  percent,  i.e.  from  8  percent  to  3  percent. 

For  purposes  of  determining  what  combinations  of  treatment 
rate  and  storage  capacity  would  meet  the  system  requirements,  the 
STORM  Quantity  Analysis  portion  was  used  to  process  twenty-two  and 
one-half  years  of  hourly  rainfall  data  recorded  at  the  U.S.  Weather 
Service  station  at  Oakland  International  Airport. 

An  initial  run  was  made  with  the  STORM  model  set  at  existing 
Special  District  No.  1  treatment  and  storage  capacities,  assumed  to  be 
0.0068  inches/hour  and  0.017  inches,  respectively,  in  excess  of  average 


Larkspur 


\  /  4  V-  >Wartme^  ’  •  * 

S/  •  xfTp^r  '  v  •  >x.  i  •  r\ 

j>J  J  Ai/Tl  ;  r  ., 

_ T  ]>  /  ^  7  ^tl.Soor^Kii®  7-  C  f  ^  /  42^^'  J"  I 

-r^  K?  '  7^.  ^\=j  1  <■  v  ^  /  v  '"*  ' s 


:.  i 


R'RiW.,, 

BrooksT^  <* 

lsl*"d  .  ;V^#.'.-  >-\-A\^' 


S*n  ^Wa  P*9 

\ 


Pfiones  b«) 


Point  Bor>,t*\[ 


eWedivT'U>®Ur°n 

%,  ^W^Un, 
^sAsalito  V_ _ Jlsl< 


mLBANY%\  .*>/*-■>■ 


Pl*  r>v 


emeryvillexX^^^x  S\v#  v 


’■**  -  U1  3 


'it  \ 

KUnd  \  t 


^PIEDMONT 


‘'"'■'VMoraga^ 


treatment  plant 


[BMUDINTERCEPTOR 


’  srv*; 


I  mtt»  Mai  ./V* 

v‘'>  Afuii**..  /.,f«U 


xf> 

San  Francisco  "yr 

Internationa!  Atrporl 


v  \  xWfo,,f5C 


I . 


v  \\v  ' 


)  , 


v.  .  ..  FIGURE  u 

Vicinity  Map  of  EBMUD  Servi 


ce  Area 


dry  weather  flow  requirements.  For  this  run,  an  existing  Di strict-wide 
gross  infiltration  rate  of  11.1  percent  of  rain,  developed  in  an  earlier 
study  [9],  was  used.  This  infiltration  rate  includes  the  effect  of  combined 
sewers  which  drain  approximately  4  percent  of  the  total  area.  The  run 
showed  an  average  incidence  of  10.  9  overflows  per  year,  which  is  in  good 
agreement  with  historical  data. 

All  subsequent  runs  of  STORM  were  based  on  the  assumption 
that  all  combined  sewers  were  separated  from  the  sanitary  system. 
T  reatment- sto  rage  combinations  were  examined  for  three  alternative 
infiltration  rates  (expressed  in  percent  of  Oakland  Airport  rainfall): 

1.  8%  -  Gross  Infiltration  Ratio  without  combined  sewers 

2.  5%  -  Removal  of  combined  sewers  plus  80%  of  "direct 

connection  s" 

3.  3%  -  Additional  reduction  by  removal  of  50%  of 

"percolation  infiltration," 

In  order  to  provide  enough  points  for  curve  plotting,  a  total 
of  49  infiltration-treatment- storage  combinations  we  re  run  as  shown  in 
Table  13.  Treatment  rates  ranging  from  0.003  to  0.03  inches  per  hour 
(107.  5  to  1075  MGD)  and  storage  capacities  ranging  from  0,001  inches 
to  0.08  inches  (1.5  to  120  million  gallons)  were  analyzed. 

Table  1  3  summarizes  the  events  per  year,  number  of  overflows 
per  year  and  quantity  of  overflow  per  year  for  each  of  the  three  assumed 
infiltration  rates. 

For  each  of  the  three  infiltration  rates,  the  average  number 
of  overflows  per  year  are  shown  graphically  in  Figure  15,  As  expected, 
the  number  of  overflows  decreases  as  the  storage  capacity  or  treatment 
rate  is  increased  or  as  infiltration  is  reduced  by  upstream  extraneous 
control  measures,  as  represented  by  the  reduction  in  infiltration  rates. 
Storage  required  to  totally  contain  the  inflow  from  the  22  year  period 
of  rainfall  record  is  shown  in  Figure  16.  It  is  not  possible  to  recommend 
the  optimum  combination  of  treatment  rate  and  storage  based  solely  on 
the  number  of  overflows.  Nor  is  this  parameter  sufficient  for  use 
in  selecting  storms  for  further  analysis  in  the  transport  model  phase. 

The  quantity  of  overflow  on  an  average  annual  basis  is  shown 
in  Figures  17,  18,  and  19  for  each  of  the  assumed  infiltration  rates. 

The  ave  rage  quantity  of  ove  r flow  al so  dec  rease s  as  the  number  of  overflows 
decrease.  The  relationships  between  quantity  and  number  of  overflows 
are  depicted  in  the  figures. 

Figures  17,  18,  and  19  show  the  relative  effectiveness  of 

treatment  and  storage  in  reducing  the  number  and  volume  of  overflows, 
and  will  serve  as  input  for  further  analysis  in  which  the  economics  of 
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treatment- storage  combinations  and  the  overflow  criteria  likely  to  be 
required  by  the  State  and  Regional  Water  Quality  Control  Boards  are 
conside  red. 
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STORAGE  REQUIREMENT- INCHES 


TREATMENT  RATE -MOD 
200  400  600  800  1000 


FIGURE  16 

Zero  Overflows  Requirements 


STORAGE  REQUIREMENT  -  MILLIONS  OF  GALLONS 
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AVERAGE  ANNUAL  OVERFLOW,  MILLIONS  OF  GALLONS 
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FIGURE  18 

Overflow  Volumes  at  5  Percent  Infiltration  Rate 
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AVERAGE  ANNUAL  OVERFLOW,  MILLIONS  OF  GALLONS 
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FIGURE  19 

Overflow  Volumes  at  3  Percent  Infiltration  Rate 
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AVERAGE  ANNUAL  OVERFLOW,  MILLIONS  OF  GALLONS 


SECTION  7 

APPLICATION  OF  "STORM"  -  EXAMPLES 


The  purpose  of  this  section  is  to  present  some  examples  of 
other  possible  applications  of  STORM  in  addition  to  those  described  in 
Section  6.  Four  applications  will  be  shown: 

1 .  Computation  of  the  quantity  of  storm  runoff  by 
month  and  for  single  storm  events. 

2.  Computation  of  pollutographs  for  single  storm 
events. 

3.  Use  of  STORM  to  find  the  most  economical 
treatment- storage  combinations  that  meet  system 
overflow  constraints. 

4.  Analysis  of  changes  in  the  quantity  and  quality  of 
urban  runoff  due  to  alternative  land  use 
management  schemes. 

A  prototype  area  was  selected  for  these  applications  of 
STORM,  i.e.,  the  Castro  Valley  watershed  near  Oakland,  California. 
Figure  20  shows  the  USGS  map  of  the  area  with  the  watershed  boundary 
and  the  location  of  the  watershed  relative  to  the  San  Francisco  Bay  area. 
Some  land  use  information  can  be  obtained  from  the  map  but  an  aerial 
photograph  plus  ground  reconnaissance  were  necessary  to  obtain  a  better 
understanding  of  land  uses  in  the  basin.  Table  14  gives  a  summary  of 
the  estimated  hydrogeometric  characteristics  of  the  watershed.  Runoff 
statistics  for  the  watershed  were  generated  by  processing  hourly  rainfall 
data  for  the  17-month  period  from  November  1971  through  March  1973. 


COMPUTATION  OF  THE  QUANTITY  OF  STORM  RUNOFF 

In  the  initial  application  of  STORM  to  Castro  Valley,  the 
program  was  calibrated  to  give  the  best  comparison  between  computed 
and  observed  values  of 

1.  Average  annual  precipitation; 

2.  Average  annual  runoff; 

3.  Monthly  runoff  volumes;  and 

4.  Individual  storm  event  volumes. 

It  is  not  appropriate  to  make  comparisons  with  the 
instantaneous  measurements  of  discharge  because  the  program  computes 
runoff  as  hourly  volumes  only.  The  hourly  volumes  should  reflect  the 
general  shape  of  the  observed  hydrograph  and  its  volume,  although  the 


-56- 


QU*0**ANG.  <  <  Of. A ' 


HAYWARD.  CALIF. 


observed  hydrograph  will  tend  to  lag  behind  the  computed  hydrograph  in 
larger  basins  where  the  time  of  concentration  is  greater  than  one  hour. 
The  program  can  be  applied  to  larger  basins  where  this  lag  problem  exists 
as  long  as  one  realizes  the  impact  on  the  analysis.  In  storage  analyses, 
however,  this  problem  is  generally  not  critical. 

TABLE  14 

Hydrogeometric  Data 


Land  Use 

Percent 
of  Area 

Percent 

Impervious 

Length  of 
Street  Gutters 

Single  Family  Residential 

70 

40 

275( ft/ac ) 

Multiple  Family  Residential 

3 

50 

430 

Commercial 

7 

80 

400 

Open  or  Park 

20 

2 

20 

Area  of  the  watershed  =  3136  acres 
Depression  Storage  =  .  10  inch 


Depression  Storage  Recovery,  inches/day 


January 

0.  05 

July 

0.  28 

F  ebruary 

0.  07 

August 

0.  25 

March 

0.  12 

September 

0.  20 

April 

0.  17 

October 

0.  13 

May 

0.  23 

November 

0.  07 

June 

0.  26 

December 

0.05 

Runoff  Coefficient  for  Pervious  areas  =  0.45 
Runoff  Coefficient  for  Impervious  area  =  0.90 

Area  weighted  basin  average  runoff  coefficient  =  0.  61 


The  program  parameters  which  were  "tuned"  in  the  calibration 
process  were: 

1.  Rainfall  factor  relating  basin  average  rainfall  to 
the  gage  rainfall; 

2.  Depression  storage  and  the  rate  of  recovery  of 
depression  storage; 

3.  Imperviousness  of  the  land  uses;  and 

4.  Pervious  and  impervious  area  runoff  coefficients. 
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The  recording  rain  gage  for  Castro  Valley  is  centrally  located  in  the 
watershed  and  was  assumed  to  reflect  basin  average  precipitation. 
Calibrated  values  of  the  other  parameters  are  those  listed  in  Table  14. 

Table  15  shows  the  computed  and  observed  data  for  monthly 
runoff  and  for  total  runoff  over  the  17-month  period.  Computed  average 
annual  rainfall  and  runoff  information  is  obtained  directly  from  the 
program's  average  annual  summary.  Monthly  volumes  were  computed 
from  the  EVENT  output.  Figures  21  and  22  show  the  computed  and  observed 
hydrograph  fit  for  two  events  in  the  17-month  record. 


TAB  EE  15 

Coaparlson  of  Observed  and  Computed 
Monthly  Runoff  Volunas  fro*  Castro  Valley,  California 


Year 

Month 

Observed  Runoff 
(Inches)* 

Computed 

Runoff 

(Inch**)** 

1971 

1  1 

.  35 

.73 

12 

1.  50 

2.  29 

1972 

1 

.  51 

.  58 

2 

.  60 

.  58 

3 

, 1 7*** 

0 

4 

.  34 

.44 

5 

,  1 0  * * 

0 

6 

.  17 

.  16 

7 

. l  l*** 

0 

8 

, 08**# 

0 

9 

.  21 

.48 

10 

1.01 

1.  51 

1  1 

2.  85 

2.  82 

12 

.  91 

1.03 

1973 

l 

5.78 

4.  13 

2 

4.  34 

2.45 

3 

2.01 

1 .71 

TOTAL  21.00 

18.  90 

11.70  bascf low 
TOTAi.  -  20.  60 

♦From  fJSGS  Records. 

-■-Does  not  include  baseflow.  Baseflow  is  approximately 
0.  10  inches  per  month. 

♦  ♦♦No  rainfall  recorded. 
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DISCHARGE  Cfs 


CASTRO  VALLEY 


HOUR  OF  DAY 
(10/11/72) 


FIGURE  21 

Comparison  of  Observed  Hourly  Values  of  Runoff 
With  Values  Computed  by  "STORM" 
November  11,  1972 
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CASTRO  VALLEY 


HOUR  OF  DAY 
(2/6/73) 


FIGURE  22 


Comparison  of  Observed  Hourly  Values  of  Runoff 
With  Values  Computed  by  "STORM" 
February  6,  1973 


n 
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COMPUTATION  OF  THE  QUALITY  OF  RUNOFF 


The  zero  storage,  zero  treatment  combination  was  also  used 
for  an  initial  calibration  of  the  quality  portion  of  STORM.  The  quality  cali¬ 
bration  was  made  difficult  because  of  the  lack  of  adequate  data.  It  is 
usually  not  economically  feasible  to  monitor  every  runoff  event,  thus 
monthly  or  average  annual  data  generally  does  not  exist.  In  most  instances, 
as  was  the  case  for  the  Castro  Valley  data,  the  quality  calibration  must 
be  made  on  data  from  several  individual  runoff  events.  However,  unless 
care  is  taken  to  obtain  a  good  sampling  over  the  duration  of  the  entire 
event,  the  comparison  of  total  amounts  of  pollutant  washoff  will  not  be 
possible.  This  was  also  the  case  in  Castro  Valley;  consequently,  most 
of  the  calibration  had  to  rely  on  a  few  data  points  per  event.  It  would 
be  highly  desirable  to  have  enough  measurements  during  a  runoff  event 
to  be  able  to  trace  the  hourly  performance  of  the  pollutant  washoff 
function. 


The  quality  calibration  for  Castro  Valley  was  made  on  the  basis 
of  comparisons  of  computed  and  observed  concentrations  for  individual 
events.  The  parameters  calibrated  were: 

1.  Dust  and  dirt  accumulation  rates; 

2.  Pollutant  composition  of  the  dust  and  dirt;  and 

3.  The  exponent  in  the  pollutant  washoff  equation. 

The  initial  pollutant  loading  rates  used  in  the  calibration  were 
those  from  the  Chicago  APWA  study  [5]  that  was  discussed  in  Section  3. 
These  data  are  programmed  as  default  values  to  be  used  if  other  values 
for  these  parameters  are  not  specified  as  input  data.  Table  16  shows 
the  calibrated  values  of  the  pollutant  loading  and  washoff  parameters. 
A  comparison  of  these  values  with  the  default  values  listed  in  Section  5 
reveals  that  the  Dust  and  Dirt  accumulation  rates  had  to  be  increased 
by  a  factor  of  two  except  for  the  open  or  park  area  which  was  increased 
by  a  factor  of  about  six.  Pollutant  composition  of  the  dust  and  dirt  was 
also  increased  by  a  factor  of  four  for  all  parameters  except  suspended 
and  settleable  solids. 


TABLE  16 

Pollutant  Loading  and  Washoff  Parameters 


Pounds  Pollutant 


Land  Use 

Dust  and  Dirt 
Accumulation 
lbs/day/100  ft.  gutter 

per  100  lbs  Dust  and  Dirt 
Susp.  Sett. 

Solids  Solids  BOD  N 

P04 

Single  Family 

1.4 

11.1 

1.  1 

2.  0 

.  19 

.  02 

Multiple  Family 

4.  6 

8.0 

.  8 

1.44 

.  24 

.  02 

Commercial 

6.  6 

17.  0 

1.7 

3.08 

.  16 

.  03 

Open  or  Park 

9.  2 

11.1 

1.  1 

2.0 

.  19 

.  02 

Street  Sweeping  Efficiency 

=  70% 

Washoff  exponent  =  2,0 


Comparisons  of  computed  and  measured  BOD  concentrations 
are  shown  in  Figure  23.  Agreement  in  the  first  case  is  very  good,  however 
computed  values  for  the  second  storm  are  consistently  low.  Notice  that 
the  first  storm  occurred  in  November  of  1972  which  is  the  beginning 
of  the  winter  rainy  season.  The  good  agreement  in  this  case  indicates 
that  the  accumulation  of  pollutant  loads  over  the  dry  summer  period  is 
being  computed  well  by  STORM. 

The  second  comparison  shown  was  for  a  storm  in  February 
which  is  near  the  end  of  the  winter.  Since  the  computed  values  are 
consistently  low,  the  implication  is  that  the  pollutant  load  on  the  watershed 
at  the  beginning  of  the  storm  was  too  small.  This  is  quite  likely  the  case 
because  during  the  period  between  the  November  storm  and  the  February 
storm,  the  computer  program  has  been  performing  a  constant  accounting 
of  the  pollutant  load  on  the  watershed,  i.e.,  how  much  was  there  at  the 
beginning  of  each  storm,  how  much  was  left  at  the  end  of  each  storm, 
and  how  much  accumulated  between  storms.  If  we  assume,  on  the  basis 
of  Figure  23  (the  November  storm),  that  the  pollutant  load  at  the  beginning 
of  the  rainy  season  is  correct  and  that  the  computation  of  washoff  rate 
is  correct,  then  the  rate  of  pollutant  accumulation  (dust  and  dirt  accumu¬ 
lation)  on  the  watershed  must  be  larger  for  the  prototype  than  the  rate 
being  used  in  the  model.  As  a  result,  the  watershed  is  being  washed 
overly  clean  by  STORM  during  the  rainy  season. 

The  next  step  in  the  calibration  procedure,  therefore,  would 
be  to  increase  the  daily  rate  of  dust  and  dirt  accumulation  during  the 
winter  period  which  would  result  in  larger  pollutant  loads  at  the  end  of  the 
rainy  season. 

STORAGE-TREATMENT  ANALYSIS  EXAMPLE 

The  primary  purpose  of  the  STORM  program  is  to  analyze  the 
effectiveness  of  storage  and  treatment  facilities  for  use  in  controlling  the 
quantity  and  quality  of  urban  storm  runoff.  The  criteria  for  control  of 
the  storm  runoff  may  be  in  terms  of  maximum  allowable: 

1.  Overflow  events  per  year: 

2.  Volume  of  overflow  per  year; 

3.  Volume  of  overflow  during  some  design  storm  event; 

4.  Pounds  of  BOD  (or  other  pollutant)  overflow  per 
year;  and/or 

5.  Pounds  of  BOD  (or  other  pollutant)  overflow  during 
some  design  storm  event. 

It  may  be  possible  to  achieve  the  desired  control  of  runoff 
through  a  number  of  different  combinations  of  storage  and  treatment. 

A  cost  analysis  is  then  necessary  to  determine  which  storage-treatment 
alternative  achieves  the  desired  control  at  least  cost. 


mg/I  BOD  mg/I 
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HOUR  0/  DAY 

(2/5/73)  (2/6/73) 


FIGURE  23 

Comparison  of  Observed  Hourly  values  of  BOD 
Concentration  With  Values  Computed  by  "STORM” 
November  11,  1972  and  February  6,  1973 
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The  decision  criteria  for  the  purposes  of  this  hypothetical 
example  are  as  follows: 

1.  No  more  than  five  overflows  per  year;  and 

2.  No  more  than  10,000  pounds  of  BOD  overflow 
pe  r  year. 

In  order  to  determine  what  sto  rage  -  treatment  rate  combi¬ 
nations  will  meet  these  objectives,  the  STORM  program  was  run  for 
several  different  storages  for  each  of  several  treatment  rates.  For  each 
storage-treatment  combination,  information  on  the  average  annual  number 
of  overflows,  inches  of  overflow,  and  pounds  of  BOD  (or  other  pollutant) 
was  obtained  directly  from  the  output  (see  Tables  8  and  10  for  example). 

This  information  canbe  plotted,  such  as  in  Figure  10  of  Section 
f>,  to  enhance  the  analysis.  For  this  problem,  however,  it  would  be  better 
to  first  plot  Number  of  Overflows/year  vs.  Storage  Capacity  and  Pounds 
of  BOD  Overflow/Year  vs.  Storage  Capacity,  both  plots  showing  lines  of 
equal  treatment  rates.  The  two  plots  are  shown  in  Figure  24.  Notice 
on  the  BOD  overflow  curves  that  the  overflow  frequency  is  also  shown. 
It  can  be  seen  from  this  figure  that  overflow  frequency  is  a  more  severe 
requirement  for  the  system  if  its  storage  capacity  is  less  than  0.3  inches. 
For  greater  storage  capacities,  however,  the  quality  constraint  is  more 
limiting.  Thus  all  treatment- storage  combinations  shownbelow  the  shaded 
line  in  Figure  25  are  acceptable  from  the  standpoint  of  meeting  or  exceeding 
the  system  performance  criteria.  However,  it  can  be  seen  that  for  a 
given  treatment  rate,  the  smallest  allowable  storage  will  be  that  identified 
on  the  shaded  performance  line  in  Figure  25.  Thus,  in  our  example 
the  most  economical  control  system  that  meets  the  performance  criteria 
will  be  one  of  the  following  three: 

1.  Treatment  =  0.01  in/hr  Storage  =  0.46  in 

2.  Treatment  =  0,03  in/hr  Storage  =  0.28  in 

3.  Treatment  =  0.05  in/hr  Storage  =  0.20  in 

Economic  Analysis  of  these  three  alternatives  would  identify  the  most 
economical  system. 


LAND  USE  MANAGEMENT  ANALYSIS 

Another  possible  application  of  STORM  is  to  estimate  the 
impact  of  proposed  land  use  changes  within  a  watershed  on  the  quantity 
and  quality  of  urban  runoff.  To  illustrate  this  use,  it  was  assumed  that 
the  upper  arm  of  the  Castro  Valley  wate  rshed  (see  Figure  20)  would  be 
developed.  This  470  acre  area  is  presently  open  space.  Under  the 
assumed  development  plan,  single  family  residential  dwellings  would  be 
constructed  on  67  percent  of  the  area.  The  remaining  33  percent  of 
the  area  would  be  developed  as  multiple  family  residences. 
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FIGURE  24 

Relationship  of  Overflow  Frequency  and  Pollutant 
Overflows  to  Treatment-Storage  Alternatives 
Castro  Valley 
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FIGURE  7 5 

Most  Economical  Treatment-Storage  Combinations 
To  Meet  System  Overflow  Constraints 


Table  17  shows  the  results  of  applying  STORM  to  the  470  acre 
subarea  in  both  its  present  undeveloped  state  and  in  the  proposed 
developed  state.  In  the  developed  state,  the  annual  quantity  of  runoff 
can  be  expected  to  increase  by  40  percent.  The  annual  BOD,  on  the 
other  hand,  will  increase  more  than  400  percent. 

The  impact  of  the  proposed  development  on  the  watershed  as 
a  whole  is  shown  in  Table  18.  As  expected,  this  impact  is  significantly 
less.  As  a  result  of  the  proposed  development,  annual  storm  runoff  from 
the  watershed  can  be  expected  to  increase  by  approximately  5  percent. 
The  annual  BOD  load  washed  off  the  watershed  will  increase  approxi¬ 
mately  14  percent. 


TABLE  17 

Effect  of  Changing  Land  Use  on 
Storm  Runoff  from  470  Acre  Subarea 


1 - 

Annual  Runoff 

Annual  BOD  Load 

Land  Use 

Inches 

Pounds 

Existing 

8.  34 

5,  700 

100%  open 

Proposed 

11.72 

23,  400 

67%  single  family  residential 

3  3%  multiple  family  residential 

TABLE  18 

Effect  of  Changing  Land  Use  of  470  Acre  Subarea  on  Storm  Runoff 
From  Entire  Castro  Valley  Watershed  (3140  Acres) 


Annual  Runoff 

Annual  BOD  Load 

Land  Use 

Inches 

Pounds 

Existing 

70%  single  family  residential 

3%  multiple  family  residential 

7%  commercial 

20%  open 

11.08 

130, 500 

Proposed 

80%  single  family  residential 

8%  multiple  family  residential 

7%  commercial 

5%  open 

1 1.  59 

148, 300 

SECTION  8 

TRANSFERABILITY  OF  WATER  QUALITY  DATA 


VARIED  INTERPRETATIONS 

The  most  difficult  problem  encountered  in  transferring  water 
quality  data  measured  in  one  area  to  another  area  that  has  no  data  is 
that  the  data  is  usually  presented  in  such  a  way  that  it  is  difficult  or 
impossible  to  separate  the  hydrologic  effects  from  the  pollution  data. 
As  was  pointed  out  in  Section  3,  the  common  denominator  for  trans¬ 
ferring  urban  runoff  quality  data  from  one  area  to  another  is  the  rate 
of  accumulation  of  pollutants  on  the  watershed  as  a  function  of  land  use 
and  some  length  or  area  parameter,  i.e.,  lbs/year/acre,  lbs/day/acre, 
Ibs/day/mile  of  street,  etc.  Much  of  the  reported  data,  however,  is  in 
terms  of  pounds  of  pollutant  washed  off  the  watershed/acre/inch  of  rainfall 
(which  does  not  correlate  with  inches  of  rainfall),  or  in  terms  of  mean 
concentrations  of  the  runoff  from  the  area.  If  sufficient  hydrologic  data 
is  presented,  this  data  can  be  reduced  to  its  basic  rate  of  buildup  form, 
but  it  is  a  tedious  task. 

Even  if  the  data  can  be  reduced,  it  is  likely  that  it  reflects 
cumulative  runoff  effects  from  several  different  land  use  areas.  In  such 
cases  the  loading  rates  for  each  land  use  could  be  derived  by  developing 
weighting  factors  based  on  the  total  length  of  gutters,  each  land  ■  se 
classification  and  the  APWA  data  presented  in  Section  3.  An  example 
of  this  weighting  method  is  detailed  below. 

Sometimes  yearly  stormwater  mass  emissions  are  reported, 
i.e.,  lbs/year  of  a  pollutant  discharged  from  a  watershed.  Loading  rates 
in  terms  of  land  use  can  also  be  derived  from  this  type  of  data  using 
weighting  factors. 


EXAMPLE  OF  DERIVATION  OF  POLLUTANT  ACCUMULATION 
RATES  FROM  MASS  EMISSION  DATA 

Weighting  factors  are  formed  for  each  pollutant  as  the  product 
of  the  D/D  rate  x  pollutants  per  unit  weight  of  D/D  x  total  gutter  length 
in  the  land  area  of  the  given  land  use.  (See  Table  19  for  BOD  Weighting 
F actors.  ) 
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TABLE  19 

BOD  Weighting  Factors  for  Reducing  Total  Mass  Emission 
Rates  from  a  Mixed  Land  Use  Watershed  to  Rates  of 
Accumulation  Per  Land  Use  Area 


Land  Use 

Weighting  Factor* 

Single  Family  Residential 

W1 

=  0.7  x  5,0  x  G^ 

Multiple  Family  Residential 

W2 

=  2.  3  x  3.  6  x  G2 

Commercial 

W3 

=  3.3  x  7,7  x  Gj 

Industrial 

W4 

=  4.  6  x  3. 0  x  Ga 
4 

Open  or  Park 

W5 

=  1.5  x  5.0  x  G^ 

*G  is  the  total  length  of  gutters  (1000's  of  feet)  in  the  portion  of 
the  watershed  of  the  given  land  use  type.  j 


To  illustrate  the  use  of  these  factors,  assume  that  it  is  reported 
that  the  BOD  mass  emission  rate  from  a  watershed  is  2,  500  pounds  per 
month,  and  that  the  watershed  has  the  following  pertinent  features: 


Land  Use  Area,  Acres 

Single  Family  Residential  197 

Multiple  Family  Residential  223 

Industrial  686 


Total  Length  of  Gutters, 
1000's  of  feet 

go 

100 

160 


The  BOD  weighting  factors  for  this  watershed  are: 

Wj  =  0.7  x  5.0  x  80  =  280 

W2  =  2.  3  x  3.  6  x  100  =  830 

W3  =  4.  6  x  3.  0  x  160  =  2210 

Total  3320 


For  a  30-day  month,  the  daily  BOD  accumulation  rate  by  land  use  area 
is  thus: 


Single  Family  Residential  (2500/30)  x  r  7.0 

Q  •>  rt 

Multiple  Family  Residential  (2500/30)  x  “20.8 

Industrial  (2500/30)  x  -  55.  5 


Total  =  83  lbs/day 


We  can  assume  that  the  APWA  valueB  of  the  amount  of  pollutant  per  unit 
of  D/D  apply  to  the  area.  These  values  are  then  used  to  derive  the 
D/D  rates  on  each  watershed.  However,  recall  that  the  values  for  BOD, 
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N,  and  P04  are  only  for  the  soluble  fraction  of  the  constituent.  Since  the 
model  assumes  that  non-soluble  BOD  is  10%  of  the  suspended  solids 
load  (we  neglect  the  2%  contribution  from  settleable  solids),  we  must 
subtract  these  amounts  from  the  total  BOD  loads  computed  above  before 
computing  the  D/D  rate.  If  the  suspended  solids  loads  for  the  Single 
Family.  Multiple  Family  and  Industrial  areas  are  23,  47,  and  200  lbs/day, 
respectively,  the  soluble  BOD  loading  rates  are: 

Single  Family  Residential  7.0  -  2.3  -  4.7  lbs/day 

Multiple  Family  Residential  20.  8  -  4.  7  =  16.  1  lbs/day 

Industrial  55.5  -  20.0  =  35.5  lbs/day 

Finally,  dividing  these  values  by  the  pounds  of  BOD  per  pound  of  D/D 
and  by  the  total  gutter  length  in  each  area  gives  the  unit  rates  of  D/D 
accumulation  for  the  watershed: 


Land  Use 


Single  Family  Residential 

Commercial 

Industrial 


D/D  Accumulation  Rate, 
lbs/day/ 100  ft.  of  gutter 

4. 7/0. 0050/800  =1.2 
15. 1/0.0036/1000  =4.4 
35. 5/. 0030/1600  =7.4 


CAVEAT 


A  substantial  amount  of  data  exists  on  the  quality  and  quantity 
of  urban  runoff.  References  10  and  11  include  fairly  complete  citations 
of  the  data  sources  that  exist.  As  previously  stated,  however,  much 
of  the  data  are  presented  in  a  form  that  is  inappropriate  for  translation 
to  other  areas. 

For  data  that  can  be  reduced  to  its  basic  form  (rate  of  accumu¬ 
lation  on  the  watershed  by  land  use)  two  cautions  are  in  order.  First, 
much  of  the  reported  data  is  measured  in  combined  sewer  systems.  Where 
this  is  the  case,  the  sewage  contribution  to  the  pollution  load  must  first 
be  subtracted  from  the  total  load  before  deriving  watershed  pollutant 
accumulation  rates.  Secondly,  "typical”  residential  and  industrial  areas 
in  one  locality  may  be  atypical  in  other  areas.  Before  tranferring  data 
from  these  areas  to  another  locality,  the  geometric  features  of  the  two 
areas  should  be  compared  and  (if  appropriate)  adjustments  should  be 
made  in  the  data  before  applying  it  to  a  new  area. 

Other  factors  which  undoubtedly  affect  the  rate  of  pollutant 
accumulation  are  the  general  air  pollution  characteristics  of  the  area  and 
the  general  climatology.  Pollutant  accumulation  rates  are  probably  much 
different  in  industrialized  Chicago  than  in  resort  towns.  They  would  also 
be  expected  to  differ  between  the  coastal  areas  of  California  and  the  arid 
cities  of  Arizona.  No  general  guidelines  can  be  presented  to  account 
for  such  effects. 
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